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JARED SPECK* 



Abstract. In this article, we study small perturbations of the family of 
Friedmann-Lemaitre-Robertson- Walker cosmological background solutions to 
the 1 + 3 dimensional Euler-Einstein system with a positive cosmological con- 
stant. These background solutions describe an initially uniform quiet fluid of 
positive energy density evolving in a spacetime undergoing accelerated expan- 
sion. Our nonlinear analysis shows that under the equation of state p = c^p, 
< < 1/3, the background solutions arc globally future-stable. In partic- 
ular, we prove that the perturbed spacetime solutions, which have the topo- 
logical structure [0, oo) X T'', are future causally gcodcsically complete. These 
results are extensions of previous results derived by the author in a collabo- 
ration with I. Rodnianski, in which the fluid was assumed to be irrotational. 
Our novel analysis of a fluid with non-zero vorticity is based on the use of 
suitably-deflned energy currents. 
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1. Introduction 

The Euler-Einstein system models the evolution of a dynamic spacctimcQ [Ai , 5^11/) 
containing a perfect fluid. As is discussed below, in applications to cosmology, an 
"additional term" of the form Ag^,y is often added to the equations in order to 
alter the nature of solutions. The constant A, which was first contemplated by 
Einstein |Einl7| . is known as the cosmological constant. For physical reasons to be 
discussed below, we assume throughout the article that A > is a fixed constant, 
and that the fluid equation of state is p = c^p, where p is the fluid pressure, p is the 
proper energy density, and the non-negative constant Cs is the speed of sound. As 
is fully discussed in Section [3l under these assumptions, the Euler-Einstein system 
comprises the equations^ 

(1.0.1a) Ric^,.-^i?5M'' + A.g^^=r^^, (/^, = 0, 1, 2, 3), 

(1.0.1b) -D„r"^ = 0, (Ai = 0,l,2,3), 

where Rie^j^ is the Ricci curvature tensor, R ~ ^"^RicQ^ is the scalar curvature, 
Tfii, — {p-\~p)u^u^ +pgfii' is the energy-momentum tensor of a perfect fluid, and u^^ 
is the four-velocity, a future-directed, timelike vectorfield that is subject to the unit 
normalization constraint ga/su'^u^ = — 1. The fundamental unknowns in the above 
system may be considered to be [M.,g^^,p,u^). Although we limit our discussion 
to the physically relevant case of 1 -I- 3 dimensions, we expect that our work can be 
easily generalized to apply to the case of 1 -I- n dimensions, n > 3. 

Our choice of the equation of state p = c^p, which is often made in the cosmology 
literature, is sufficient to close the above system of equations. We remark that 
without specifying an equation of state, there would be too many fluid variables, 
and not enough fluid equations. There are far too many references on this subject 
to provide a comprehensive list, but for an overview of the field, readers may consult 
e.g. the mathematical references |Car01| . jRenOSj . |Wal84j . and the cosmological 
references |PR03| , |Sah04| . As is explained in Section |4l under the equation of 
state p — c^p, there exists a well-known family of Friedmann-Lemaitre- Robertson- 
Walker (FLRW) solutions to (|1.0.1a|) - (jl.0.1b|) that are frequently used to model a 
fluid-filled universe undergoing accelerated expansion; these are the solutions that 
we investigate in detail in this article. We remark that the accelerated nature 
of the expansion in the FLRW family is generated by the positivity of 
A; when A = 0, the FLRW expansion is no longer accelerated. 

The cases p = and p = (l/3)p, which are known as the "pressureless dust" 
and "radiative" equations of state, are of special significance in the cosmological 
literature. The latter is often used as a simple model for a "radiation-dominated" 
universe, while the former for a "matter-dominated" universe. Unfortunately, as 
we will see, these two equations of state lie just outside of the scope of our main 
results, which are restricted to the parameter range < Cg < -\/l/3. We expect 
that an analogous stability result can be proved in the case Cs = 0, i.e., that the 



A spacetime is defined to be a 4— dimensional time-orientablc Lorcntzian manifold Ai together 
with a metric g^i/ on M of signature (— , -|-, +, +). 

"^Throughout the article, we work in units with SttG = c = 1, where c is the speed of light 
propagation in Maxwell's theory of electromagnctism, and G is Newton's universal gravitational 
constant. 
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failure of our proof in the case of the pressureless dust is only due to the fact that 
our methods degenerate: in the case = 0, it can be checked that the energy- 
density loses one degree of differentiability, which requires a re- working of all of the 
estimates. We will treat this case in detail in an upcoming article. In contrast, 
there may be instability in the cases Cg > -^/l/S. This is an interesting possibility 
worthy of further investigation. We state our main results roughly here; they are 
are stated more precisely as theorems in Sections [11] and [T^l 

Main Results. If < Cg < then the FLRW background solu- 

tion ([0,oo) X T^,gf,^,p,u^'), in,i^ = 0,1,2,3), to (jl.O.laj) - (jl.O.lbp . which 
describes an initially uniform quiet fluid of constant positive pressure p, 
is globally future-stable under small perturbations. In particular, small 
perturbations of the initial data corresponding to the background solution 
have maximal globally hyperbolic developments that are future^ causally 
geodesically complete. Here, g ~ — dt-^ + e^^'*' X]i=i (^^2;')^, p = g-(3+x)np^ 
and = (1,0,0,0), where >i = Sc^ and n{t) - i^/A7^)t is defined in 
()4.0.16p . Furthermore, in the wave coordinate system introduced in Sec- 
tion lS.ll the components g^i, of the perturbed metric, its inverse g'^", p, u^, 
and various coordinate derivatives of these quantities converge as t — >■ oo. 

Remark 1.0.1. Note that our results do not address the question of whether or 
not the perturbed solution is converging back to exactly the FLRW solution. 

Remark 1.0.2. In this article, we restrict our attention to spacetimes that have 
spatial slices that are diffeomorphic to T"^ (i.e., [— 7r,7r]^ with the ends identified). 
However, we expect that the local patching arguments developed in jRinOSj could 
be used to prove future-stability results for a larger class of background solutions 
featuring alternative compact spatial topologies. 

The Main Results above arc an extension of recent work by Rodnianski and 
the author [RS09j , which provides a proof of analogous results in the case that the 
fluid is irrotational. In a future article, we plan to further extend these results by 
analyzing the behavior of the fluid when the spacetime expansion is sub-exponential. 
In the case of an irrotational fluid, the fluid equations (that is, the relativistic Euler 
equations (jl.0.1b|) . which are more fully discussed in Section [3]) are equivalent to 
the following scalar equation for the derivatives 9$ of the fluid potential $ : 

(1.0.2) D„(ct"D"$)=0, 

where a = —g°'^ {da^){dp^) is the enthalpy per particle, and s = (1 — c^)/(2Cs). 
We remark that in the irrotational case, the pressure, proper energy density, and 
four- velocity can be expressed in terms of 9$ as p = -^^<^^^^ ^ P = y_^j^ o-^ and 
it^ = —a^^/^dfj_^. Before discussing the differences between our analysis of the full 
Euler system ([l.O.lbp and the scalar equation (|1.0.2I) . we provide some additional 
physical context for the problem at hand. 

Our principal motivation for inserting a positive cosmological constant into the 
Einstein equations is that there is experimental evidence for accelerated expansion. 
The first such evidence came from experiments carried out in the 1990's, which 
involved measurements of the Cosmic Microwave Background and the red shift of 



'Throughout this article, dt is future-directed. 
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type lA supernovae. The notion of accelerated expansion is perhaps best ihustrated 
through a specific example, such as the following solution to the Einstein- vacuum 
equations with a positive cosmological constant (i.e., (jl.O.f ap with T^^ = 0): the 
metric g = —dt^ + e^^* on the manifold A4 = (— cxd,cxd) x M^, where 

H = ^y A/3. In this spacetimc, the curves 7'^(t) — (t, &, 0,0), where 6 is a con- 
stant, arc (non-affinc parameterized) geodesies representing the trajectories of a 
family of physical observers. Let us denote such a curve (observer) by 7(f,). Here, 
{x^ ,x'^ ,x^) is a standard coordinate system on (—00,00) x R-^. At r = 0, the 
observers 7(0) and 7(f,) are initially separated by a distance 6, while at later times 
T, their distance is be^'^. It is roughly this kind of growth in distance at an in- 
creasing rate that leads to the term "accelerated expansion." Note that the FLRW 
background metrics from the Main Results feature this kind of behavior. Also 
note that the aforementioned distance is measured in the Riemannian manifold 
(S.r,5^^j), where Tij- = {x G A4 \ x^ = t}, and g^^^ is the Riemannian metric on 
T,r induced by _g, i.e., the first fundamental form of Et-. Consequently, it is not a 
notion inherent to spacetime itself, but instead depends on the particular foliation 

T }r G ( — cxD,oo) • 

The Main Results stated above have made implicit use of the following funda- 
mental facts concerning the Einstein equations: 

(1) The question of the local existence of solutions can be formulated as an 
initial value problem, in which suitably specified initial data give rise to 
unique local solutions. 

(2) Each initial data set satisfying the Einstein constraint equations (which are 
discussed below) launches a unique maximal solution known as the maximal 
globally hyperbolic development of the data. We emphasize that we use the 
term "maximal" to mean the largest possible globally hyperbolic0 spacetime 
determined uniquely by the data. 

These two facts, which were established by Choquet-Bruhat |CB52j and Choquet- 
Bruhat/Geroch |CBG69| respectively, are by now well-known. However, their va- 
lidity is disguised by the diffeomorphism invariance of the equations. Roughly 
speaking, this means that each spacetime "solution" is represented by an infinite 
number of solutions to the equations (|1.0.1ap - (jl.O.lbp . all of which are related 
by changes of coordinates. A closely related difficulty is that there is no canon- 
ical coordinate system for analyzing the Einstein equations. In her proof of (1), 
Choquet-Bruhat overcame these difficulties through the use of wave coordinates^ a 
special class of coordinate systems that date back at least to 1921, where they are 
featured in the work of de Bonder |dD21j . More specifically, the classic wave coor- 
dinate condition is F'' = 0, (/x = 0,1,2,3), where the F'' = g'^^T^^ (see (|3.1.2dp ) 
are the contracted Christoffel symbols of the metric. Although it is not difficult (see 
Section 1^75)) to construct a spacetimc coordinate system in which the F^ = along 
the "initial" Cauchy hypersurface E (upon which the initial data are specified), 
it is not immediately clear whether or not this condition can always be extended 
(for solutions to the Einstein equations) to include a spacetime neighborhood of E. 
Choquet-Bruhat was the first to answer this question in the affirmative. 



Recall that a spacetimc is globally hyperbolic if and only if it contains a Cauchy hypersurface 
(see Definition [322J. 

'"'A wave coordinate system is also known as harmonic gauge or de Bonder gauge. 



4 



Stability of the FLRW Family 



For the purposes of proving a local existence result, the key point is that in 
a wave coordinate system, the "gravitational part" of the Einstein equations (i.e. 
equations ()1.0.1a|) ) is equivalent to a modified system comprising quasilinear wave 
equations. Roughly speaking, the modified system is obtained by setting F*^ = 
for certain terms in (jl.O.f ap . When the modified equations are coupled to the first 
order fluid equations (jl.0.1b|) (which are well-known to be hyperbolic), the result 
is a hyperbolic system of mixed order. For such hyperbolic systems, local existence 
in a suitable Sobolcv space follows from rather standard methods based on energy 
estimates; sec e.g. [Hor97| . }SS98j . |Sog08[ , |Tay97| , and the discussion in Section 
15.51 From this perspective, the main contribution of Choquet-Bruhat's work ^CB52] 
is as follows: under suitable assumptions, the F'' = condition is propagated by 
the flow of the modified coupled system; see Section for more details. 

Before commenting further on the subtleties of the analysis, let us discuss some 
standard facts concerning the initial data. The initial data for the Euler-Einstein 
system consist of a 3— dimensional Ricmannian manifold S, and the following ten- 
sorficlds on E : a Ricmannian metric g_.^ , a symmetric two-tensor K_ji^ , a function 

p, and a vectorfield vP , (j, A; — 1,2,3). A solution consists of a 4— dimensional man- 
ifold M, a Lorentzian metric g^^, a function p, a future-directed, unit-normalized 
vectorfield u^, il^iV = 0,1,2,3), on M. satisfying (jl.O.lap - (jl.0.1b|) . and an em- 
bedding Yi ^ M such that g is the first fundamental fornjl of E, A] is the 

second fundamental forrrQ of E, the restriction of p to E is p, and (u^, u^, w'^) is the 
g— orthogonal projection of (u", u^jU^, u^) onto E; see Section [2731 for a summary of 
the conventions we use for identifying tensors inherent to E with spacetimc tensors. 
It is important to note that the initial value problem is overdetermined, and that 
the data are subject to the Gauss and Codazzi constraints: 

(1.0.3a) k - K,,t'' + it'^ab? = 2rools, 

(l-0.3b) WKa^ - t^j^b = To,\t, 

where R is the scalar curvature of g, D_ is the Levi Civita connection corresponding 
to ij. In the above formulas, indices are lowered and raised with g and g"^, and the 
future-directed unit normal TV to E is assumed to have the components iV^ = 5q . 
We note that the analysis of Section 13.2.31 implies that under the equation of state 
p = we have that Tools = ^-^Pi^oY - Pi and To^lg = ^r^puolLj, where 

uo^-^Jl + g^^uu- 

Remark 1.0.3. In this article, we do not address the issue of solving the constraint 
equations for the system (|1.0.1ap - (|1.0.1b[) . 

The fact (2) from above was settled by Choquet-Bruhat and Geroch in 1969 
jCBG69j , who showed that every initial data set satisfying the constraints (|1.0.3ap 
- (jl.0.3bp launches a unique maximal globally hyperbolic development. Roughly 
speaking, this is the largest spacetime solution to the Einstein equations that is 

^Recall that g is defined at the point x by g(X, Y) = g(X, Y) for all X,Y e Txt. 
■^Recall that K is defined at the point x by k{X, Y) =' g{DxN, Y) for all X,Y G T^t, where 
TV is the future-directed unit normal to E at x. 
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uniquely determined by the data (see Section for additional details). This is 
still a local existence result in the sense that the resulting spacetime may contain 
singularities. In particular, causal geodesies may terminate, which in physics ter- 
minology means that an observer (light ray in the case of null geodesies) may run 
into the end of spacetime in finite affine parameter. Such a singularity corresponds 
to the breakdown of the deterministic nature of the laws of physics. For space- 
times launched by initial data near that of the FLRW background solutions, our 
main result rules out the possibility of these singularities for observers (light rays) 
traveling in the "future direction." 

1.1. Comparison v^rith previous work. A noteworthy feature of our main re- 
sult is that it shows that fluids behave very differently on exponentially expanding 
backgrounds than they do in flat spacetime. More speciflcally, if one fixes a back- 
ground metric on [0, oo) x T'^ near ^^i,, then our proof can be easily adapted to show 
that the relativistic Euler equations (jl.0.1b|) on this expanding background with 
< Cs < \/l/3 have global solutions arising from data that are close to that of an 
initial uniform quiet fluid state. In contrast, Christodoulou's monograph |Chr07b| 
shows that on the Minkowski space background, shock singularities can form in 
solutions to the Euler equations arising from data that are arbitrarily close to that 
of a uniform quiet fluid state, even if the solution is irrotational. Thus, we state 
the following as a corollary of our proof: exponentially expanding spacetimes can 
stabilize perfect fluids. 

The first author to obtain global stability results for the Einstein equations with 
A > was Helmut Friedrich, first for 1 + 3— dimensional vacuum spacetimes |Fri86j . 
then later for the Einstein-Maxwell and Einstein- Yang-Mills systems [Fri91| . He 
developed a technique known as the conformal method, which translates (via a 
suitable change of variables) the question of the global stability of a solution to 
the Einstein equations into the question of local stability for the conformal field 
equations. Friedrich's work on vacuum spacetimes was later extended by Anderson 
[AndOSj to apply to 1 -t- n dimensional manifolds, where n is odd. Unfortunately, 
not all matter models seem to be amenable to the conformal method. 

An alternative to the conformal method, which was provided by Ringstrom in 
[RinOSj ■ inspired our work on the irrotational case jRS09j . In |Rin08j . Ringstrom 
used a wave coordinate approach to show the existence of an open family of future 
causally geodesically complete solutions the Einstein-scalar field system. In |Rin08] . 
the scalar field is postulated to satisfy the equation g"^ DaD^^ = ^(*&) with associ- 
ated energy-momentum tensor T^^ = (9^<I>)(9,y<i>) - [^g"^ {da^){dp^) + V{^)]gfj_^, 
where D is the Levi-Civita connection corresponding to g, and V'{^) is a non- 
linearity satisfying V{0) > 0, V'{0) = 0, V"{0) > 0. Although the cosmological 
constant is set equal to in jRin08| . these three conditions allow the nonlinearity 
V to emulate the presence of a positive cosmological constant, if $ is sufficiently 
small. In fact, the case $ = is equivalent to the vacuum Einstein equations with 
a cosmological constant A = V^(0). One of the main advantages of Ringstrom's 
framework is that it can be adapted to handle a wide variety of matter models. 
Its robustness is what prompted us to adapt it to the problem of interest in this 
article. 

Our results were further motivated by |BRR94| . in which Brauer, Rendall, and 
Reula showed a Newtonian analogue of our main result. More specifically, they 
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studied Newtonian cosmological modelf0 with a positive cosmological constant and 
with perfect fluid sources under the equation of state p = Cp'^ , where p > is the 
density, C > is a constant, and 7 > 1 is a constant. They showed that sniaU 
perturbations of a uniform quiet fluid state of constant positive density lead to a 
global solution. The fact that their work did not require the fluid to be irrotational 
strongly suggested that our results from |RS09| could be generalized to the Main 
Results stated above. 

Finally, we remark that the mechanism of stability for the problem at hand is of 
a very different nature than that encountered in the well-known body of work on the 
stability of the Minkowski spacetime solution (the original proof is }CK93j . while 
alternate proofs/extensions can be found in e.g. |BZ09| . jKN03| . |Loi09| . |LR10j . 
|SpelO| ), which features A = 0. In the latter case, there is a delicate competition 
between the size of the nonlinear terms and the dispersive nature of solutions to 
the corresponding linearized equations, while in the present case, we are aided by 
the energy-dissipative effect induced by A > 0. See |Rin08| for a more thorough 
comparison. 

1.2. Comments on the analysis. Wc now discuss our strategy for proving global 
existence. Our analysis of the metric components 5^1/ is essentially the same as in 
jRinOSj and jRS09| , but our analysis of the fluid variables differs markedly from the 
analysis of the scalar fluid equation in jRS09) , and from the analysis encountered in 
[RinOSj . Let us first discuss a simple model problem that captures the spirit of our 
analysis of the 3^,^. Consider the inhomogeneous wave equation g"^ DaDpv = F 
for the model metric g — —dt^ + e?*" Yl^a=ii'^^'^)'^ '^'^ manifold-with-boundary 
M = [0,cxd) X T'"*. Here, we are using a standard local coordinate system x^,x^,x^ 
on T''. Simple calculations imply that relative to this coordinate system, the wave 
equation can be expressed as follows: 

(1.2.1) -d'^v + e-'^*5''^dadbv = iidtvf + F. 

To estimate solutions to (|1.2.1|) . one can define the "usual" energy E{t) > by 
E-^{t) = \ Jj.sldtv)'^ + e~'^*6'^^{dav){dt,v) d'^x, and a standard integration by parts 
argument leads to the inequality 

(1.2.2) ^^E<-E+\\FU2. 

It is clear from (|1.2.2p that sufficient estimates of ||i^||L2 in terms of E will lead to 
energy decajQ. 

In Section r5.41 we decompose our modified version of the Euler-Einstein system in 
order to better see the structure of its terms. The modified system is constructed 
in a manner such that the g^i/ satisfy equations with "principal terms" that are 
similar to the terms model equation (|1.2.ip . By "principal," we are not referring to 
the degree of differentiability, but rather to the effect that these terms have on the 
solutions. As is typical in the theory of nonlinear hyperbolic PDEs, most of our 
effort goes towards analyzing the remaining "error terms," which are analogous to 
the inhomogeneous term F in the model equation. This analysis, which is carried 

^Their models were based on Newton-Cartan theory, which is a slight generalization of ordinary 
Newtonian gravitational theory that can be endowed with a highly geometric interpretation. 
^In our work below, we work with rescaled energies that are approximately constant in time. 
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out in Section 1^1 justifies the division into principal and error terms. We remark 
that the main tools used for estimating the inhomogencous terms are Sobolev-Moser 
type estimates, which we have placed in the Appendix for convenience. 

In [RS09| , our analysis of the irrotational fluid equation (|1.0.2p , was in the spirit 
of the above discussion. The new aspect compared to Ringstrom's work |Rin08| 
was that the effective metric for (|1.0.2p is not the inverse spacetime metric g^'^ , but 
is instead the reciprocal acoustical metric {m~^)^^ , an inverse Lorcntzian metric 
that features a family of sound cones as acoustical null hypersurfaces. The main 
challenge in |RS09| was that the acoustical metric can lead to instabilities if the 
pressure decays too quickly. In particular, it can degenerate if the pressure be- 
comes at a point. Addressing this difficulty was unavoidable, since the FLRW 
background solutions feature pressures that exponentially decay towards in time. 
To close our global existence argument, we had to ensure that the pressure of a 
perturbed solution decayed at the same rate as the background pressure, and in 
particular, that the pressure never becomes in finite time. 

Our analysis of the relativistic Euler equations (jl.0.1b|) is based on energy es- 
timates derived from energy currents. The energy currents (see (|6.3.13|) ) are 
vectorfields that can be used via the divergence theorem to control the evolution of 
the norms of the fluid variables and their derivatives. This approach to analyzing 
the Euler equations differs from the well-known symmetric hyperbolic framework 
(see e.g. jCII89| . |DaflOj . |Fri54| ) in that it allows the analysis to take place directly 
in the Eulerian variables (p, u) rather than an artificially introduced collection of 
unphysical state-space variables; see the introduction of |Spe09a| for a more detailed 
discussion on the merits of the energy current framework. Our energy currents ul- 
timately owe their existence to the fact that the relativistic Euler equations are a 
hyperbolic system of PDEs derivable from a Lagrangian. That is, the relativistic 
Euler equations are the Euler-Lagrange equations corresponding to a Lagrangian. 
For such hyperbolic systems, Christodoulou has developed a geometric-analytic 
framework |ChrOO| for deriving energy estimates based on the availability of energy 
currents. In fact, the use of this technique in the context of deriving energy esti- 
mates for the relativistic Euler equations was first introduced in |ChrOO| . However, 
there is a subtlety: the Lagrangian formulation of relativistic fluid dynamics de- 
scribed in [ChrOOj naturally has as its unknowns the so-called Lagrangian variables, 
while in this article, we takes our unknowns to be the Eulerian variables. In order 
to derive energy currents in the Eulerian variables, the natural way to proceed is 
to first derive them in the Lagrangian variables, and then translate them into ex- 
pressions involving the Eulerian variables. This is a decidedly nontrivial task since 
the transformation from Lagrangian to Eulerian variables is nonlinear. Further 
complicating matters is the fact that the relativistic Euler equations in Lagrangian 
variables suffer from a degeneracy that renders them just outside of the scope of 
the framework of |ChrOO| . In particular, the energy currents derived in |ChrOO| 
are not fully coercive, but are instead only semi-coercive. Nonetheless, one can 
still derive energy currents that are L^— coercive in the Eulerian variables. These 
issues are themselves of interest, and we will provide a detailed discussion of them 
in a future article. In the present article, we restrict our attention to adapting the 
Eulerian- variable energy currents used by Christodoulou in [Chr07b| . and by the 
author in |Spe09a| , |Spe09b| , without fully explaining their origin or motivating 
their cocrcivcness properties. 
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Finally, we would like to make some comments on the assumption < Cg < -^/l/S 
stated in the Main Results above. In order to do this, we have to provide 
some details concerning the behavior of the fluid norms UN-iit), Sp-p^u'-.Nit), 
and the total solution norm Sjv(Oi which are defined in (|6.I.2ap . (j6.1.2b|) . and 
(j6.1.3|) respectively. Here, > 3 is an integer. Now \i x < 1 (i.e., Cs < i/l/S), 
then the bootstrap assumption UN-iit) < e plus Sobolcv embedding together 
imply that ||u-'||ioo < Ce^^^+''-'^*e, where the small positive constant q is de- 
fined in (|8.I.2p . From this decay estimate, it follows that the length of the spa- 
tial part of w, as measured by the background metric is exponentially decay- 
ing toward the quiet background fluid state: from standard Sobolcv-Moscr es- 
timates (see Appendix |^ and the bootstrap assumption S'p_p.„»;Ar(t) < e, it 
follows that WgabU^'u^H^ = \\e^"Sabu''u''\\HN < Ce^"'E''j=^\W'\\L^h'\\H« < 
Ce^'^^* Sp-p^u';N{t) < Cee^*^*; this last estimate shows that the four-velocity 
exponentially decays towards the quiet background state, which is = 0, [j ~ 
1,2,3). In contrast, the bootstrap assumption Sp-p,u*-N{t) < e would by itself 
only lead to the estimate ||u^||icxj < Cee~^*, which would lead to the estimate 
||gabM°M^||ff« < Ce. Note that this weaker estimate does not imply that the four- 
velocity is decaying towards the uniform quiet background state. In fact, the sole 
reason we introduce the quantity UN-i(t) is to capture the additional L°° exponen- 
tial decay of the four- velocity. As we will see in Proposition 19 . 2 . Il the rapid decay 
of the fluid towards the quiet state plays a crucial role in many of our estimates. 

It is important to note that the above reasoning breaks down when x > 1 (the 
degeneracy that occurs in the case Cg = has been discussed above). We first re- 
mark that the influence of the positive cosmological constant on the behavior of the 
fluid four-velocity is essentially captured by the dissipative — 2)uju^ term on the 
right-hand side of equation (|5.7.5cp . When k > 1, equation ()5.7.5c|) suggests (we 
are assuming that A ^ is an error term) that we would at best be able to prove an 
estimate of the form ||u-'(t) H^oc < Cee~(^~^)^*; i.e., we no longer expect to have 
an improved estimate of the form < Cee^*^^"'"'^-'^* (for some q > 0), and this 

lack of improvement reverberates throughout the remaining estimates. For exam- 
ple, the gabu°'u^ term in (|9.2.14p is no longer expected to decay, so that inequality 
(|9.2.2bp must be downgraded to at best — l||^]v < CSjv- By Sobolcv embedding, 
we conclude that at best the downgraded estimate \\vP — < Ce holds. This 

weakened estimate returns to haunt us in equation (|11.1.5p . which involves the term 
2(i-i-cJJ'(m (g^g^j^)^°{tfe. This now-fatal term leads to the following downgraded 

version of inequality (lll.l.lbl) : El_^^^, .j,{t) < £;|,_p- „..jv(ii) + C E%{t) dr, in 
which we have lost the availability of the e"*^"^ damping factor. Here, Ep^p^u'-.N 
is the fluid energy defined in (|6.3.15p . and Ejv is the total solution energy defined 
in (|6.4.ip . By Proposition 110.0.21 Ep-p^^'-N ~ Sp-p^u'-N and Ejv ~ Sat for per- 
turbations of the FLRW solutions. This resulting weakened estimate allows for the 
possibility that Ep-p^u*:N{t) may grow unabatedly in time, which would eventually 
destroy the bootstrap assumption Sp-p^u'-.Nit) < e. We remark that these difficul- 
ties have nothing to do with the fact that we are studying the coupled system: they 
would remain even if we were studying the relativistic Euler equations on the fixed 
spacetime background ((— oo,oo) x T'^,^^,^). 

1.3. Outline of the structure of the paper. 
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• In Section [51 we describe our conventions for indices and introduce some 
notation for differential operators and Sobolev norms. 

• In Section [3l we introduce the Euler-Einstein system. 

• In Section |4l we use a standard ODE ansatz to derive a well-known family 
of background Friedmann-Lemaitre-Robertson- Walker (FLRW) solutions 
to the Eulcr-Einstcin system. 

• In Section [5] we introduce our version of wave coordinates and use algebraic 
identities valid in such a coordinate system to construct a modified version 
of the Euler-Einstein system. We then discuss how to construct data for 
the modified system from data for the unmodified system in a manner 
that is compatible with the wave coordinate condition. Next, wc discuss 
classical local existence for the modified system, including the continuation 
principle that is used in Section [11] We also sketch a proof of the fact 
that the modified system is equivalent to the unmodified system in a wave 
coordinate system. 

• In Section [6l we introduce the relevant norms and the related energies for 
the modified system that we use in our global existence argument. In order 
to construct the fluid energies, wc introduce the fluid equations of variation 
and the fluid energy currents. We also provide a preliminary analysis of the 
derivatives of the energies, but the inhomogeneous terms are not estimated 
until Section [9l 

• In Section [71 we introduce some bootstrap assumptions on the spacetime 
metric g^i^. We then use these assumptions to provide some linear- algebraic 
lemmas that are useful for analyzing g^i^ and the inverse metric g'^'^ . 

• In Section[8l we introduce our main bootstrap assumption, which is a small- 
ness condition on Sat, a norm of difference between the perturbed solution 
and the background solution. We also define the positive constants q and 
■Hmin, which play a fundamental role in the technical estimates of the fol- 
lowing sections. 

• Section [9] contains most of the technical estimates. We assume the boot- 
strap assumptions from the previous sections and use them to deduce es- 
timates for (7pjy, g^"^ , p, , and for the nonlinearities appearing in the 
modified equations. 

• In Section llOi we show that the Sobolev norms and energies defined in 
Section [6] are equivalent. 

• In Section [TTl we use the estimates from the previous sections to prove our 
main theorem, which is a small-data global (where "small" means close to 
the background solution) existence result for the modified equations. Our 
theorem shows that initial data satisfying the Euler-Einstein constraints, 
the wave coordinate condition, and the smallness condition lead to a future- 
geodesically complete solution of the unmodified Euler-Einstein system. 

• In Section 1121 we provide a theorem concerning the convergence of the 
global solutions as t ^> oo. 

2. Notation 

In this section, we bricfiy introduce some notation that we use in this article. 

2.1. Index conventions. Greek indices a, /?, • ■ • take on the values 0, 1, 2, 3, while 
Latin indices a, &, ■ • ■ (which we sometimes call "spatial indices") take on the values 
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1, 2, 3. Pairs of repeated indices, with one raised and one lowered, are summed (from 

to 3 if they are Greek, and from 1 to 3 if they are Latin). We raise and lower 

indices with the spacetime metric 17^^ and its inverse g^'^ . Exceptions to this rule 

include the constraint equations ()1.0.3ap - (jl.0.3bp and p.2.8a|) - ()3.2.8bp , in which 

we use the 3— metric g and its inverse g^'^ to lower and raise indices, and in 

—1^ — 

Section [121 in which all indices are lowered and raised with g^i, and g^'^ except for 

the 3— metric 3^^°'' , which has ff(^) as its corresponding inverse metric. 

2.2. Coordinate systems and differential operators. Throughout this article, 
we work in a standard local coordinate system x^,x'^,x^ on T'^. Although strictly 
speaking this coordinate system is not globally well-defined, the vectorfields dj == 
g|j are globally well-defined. This coordinate system extends to a local coordinate 
system x^,x^,x'^,x^ on manifolds with boundary of the form Ai = [0,r) x T'^, 
and wc often write t instead of In this local coordinate system, the background 
FLRW metric g is of the form (|4.0.10|) . We write to denote the coordinate 
derivative gfjr, and we often write dt instead of do- Throughout the article, we will 
perform all of our computations with respect to the fixed frame \dn] _„ , o 

If (3 = (rii, 71.2, "-3) is a triplet of non- negative integers, then we define the spatial 
multi-index coordinate differential operator ds by ds d^^ d^" . We denote the 
order of a by |q?|, where |d| ni + n2 + n^. 

We write 



(2.2.1) 

r s 
a—1 a—1 

(where T is defined in (j3.1.2dp ) to denote the components of the covariant de- 
rivative of a tensorficld on Ai with components (relative to the coordinate frame 
introduced above) T^J.'.'^^ ■ 

We write c^'-^-'T'^j.'.'.'^^ to denote the array containing of all of the TV*'' order space- 
time coordinate derivatives (including time derivatives) of the component T'^J.'.'.'^^ ■ 
Similarly, we write 

d^^^T^^l'-.Z *° denote the array containing of all A^*'' order spa- 
tial coordinate derivatives of the component T^^W.'^^ . We omit the superscript 
when = 1. 

2.3. Identification of spacetime tensors and spatial tensors. We will often 
view T'^ as an embedded submanifold of the spacetime M under an embedding tt 
of the form : ^ {t} xT^ CL M, Lt{x^ , x"^ , x^) = (t, x\ x^, x^). Note that the 
embedding is a diffeomorphism between T'^ and {t} x T'^. We will often suppress 
the embedding by identifying T'^ with its image tt(T'^). Furthermore, if T^j^.'.'.j'' is 
a T"^— inherent "spatial" tensorficld, then there is a unique "spacetime" tensorfield 
defined along ^((T^) ~ such that l^T = T and such that T' is tangent 
to tt(T'^). Here il denotes the puUback by it- We will often identify T with T', and 
use the same symbol to denote both. e.g. T^jl[[.jl — T^l'.'-'-'Z_^\ we especially apply 

"'^^Recall that T^^^. '. 'i^^ is tangent to tt(T^) if any contraction of any upstairs (downstairs) index 
with the unit normal covector (unit normal vector n'') results in 0; for downstairs indices, this 
notion depends on the spacetime metric gfiv 
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this convention to the initial data. All of these standard identifications should be 
clear in context. 



2.4. Norms. All of the Sobolev norms we use are defined relative to the local 
coordinate system cc^, x^, on T'^ introduced above. We remark that our norms are 
not coordinate invariant quantities, since we work with the norms of the components 
of tensorfields relative to this coordinate system. If / is a function defined on the 
hypersurface {x ^ AA \ t ^ const} ~ T"^, then relative to this coordinate system, 
we define the standard Sobolev norm ||/||rrjv follows, where (Px dx^dx'^dx^ : 



(2.4.1) ||/L« = fE / \daf{t,x\x^x')\'d'x 

The symbol d^x represents a slight abuse of notatior[3 since the coordinate system 
x^,x'^,x^ is not globally well-defined on T'^ (even though the differential operators 
ds are). 

Using the above notation, we can write the A^*'' order homogeneous Sobolev 
norm of / as 




(2.4.2) lla^^VlL. = E ll^-s/IL- 

If .8 C M" or C T", then Cf (j?) denotes the set of A-times continuously 
differentiable functions (either scalar or array-valued, depending on context) on 
the interior of A with bounded derivatives up to order N that extend continuously 
to the closure of A. The norm of a function F £ C^{K) is defined by 

(2.4.3) \F\N,ji= V esssup|9,-i^(.)|, 

\I\<N 

where df is a multi-indexed operator representing repeated partial differentiation 
with respect to the arguments • of F, which may be either spacetime coordinates 
or metric/fiuid components depending on context. When A^ = 0, we also use the 
notation 



(2.4.4) |FU = esssup|F(.)|. 
Furthermore, wc use the notation 

(2.4.5) |fWu= J2 

\I\=N 

The quantity is a measure of the size of the A*'' order derivatives of F on 

the set ^. In the case that = T'^ , we sometimes use the more familiar notation 



^ ^ A precise definition of the norm 1)2.4.111 would involve the use of an atlas on and a partition 
of unity that is subordinate to the corresponding covering. 
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(2.4.6) 



dof 



ess sup 



(2.4.7) 




dcf 



'N — 



\S\<N 



If A is an m X n array-valued function with entries Ajk, (1 < J < ti, 1 < fc < n), 
then in Section [9l we write e.g. ||yl||^jv to denote the m x n array whose entries 
are || A^a; || . We use similar notation for other norms of A. 

If / C M is an interval and X is a normed function space, then we use the 
notation C^{I,X) to denote the set of iV-times continuously differentiable maps 
from / into X. 

2.5. Running constants. We use C to denote a running constant that is free to 
vary from line to line. In general, it can depend on N (see (|8.0.6p ) and A, but 
can be chosen to be independent of all functions {g^^, P,u^), {jJL^v = 0,1,2,3), 
(j = 1,2,3), that are sufficiently close to the background solution {gf^„,P,u^) of 
Section 131 We sometimes use notation such as C{N) to indicate the dependence of 
C on quantities that are peripheral to the main argument. Occasionally, we use c, 
C*, etc., to denote a constant that plays a distinguished role in the analysis. We 
remark that many of the constants blow-up as A — s- 0^. 

2.6. A warning on the sign of fig. Although we often choose notation that 
agrees with the notation used by Ringstrom in |Rin08j . our reduced wave operator 
fig == g^^cladp has the opposite sign of the one in |Rin08| . It has the same sign as 
the one used in |RS09| . 



In this section, we provide a detailed introduction to the Euler-Einstcin system, 
including the notion of a perfect fluid. We then discuss several aspects of the initial 
value problem formulation of the Eulcr-Einstein system, including the initial data 
and the maximal globally hyperbolic development of the data. 

3.1. Introduction. The Einstein equations connect the Einstein tensor Ric^^ — 
■^gf_,i,R, which contains information about the curvature of the spacetime (A^, g^i/) 
to the energy-momentum tensor T^y, which models the matter content of spacetime: 

(3.1.1) Ric^^ - ^i?,g^,. + Kg^^ = T^^, 

where Ric^,y is the Ricci curvature tensor, R is the scalar curvature, and A is the 
cosmological constant. We remark that the stability results proved in this 
article heavily depend upon the assumption A > 0. Recall that the Ricci 
curvature tensor and scalar curvature are defined in terms of the Riemann curvature 
tenso^^ Riem^,^j,^. which can be expressed in terms of the Christoffel symbols F^"^ 



^■^Undcr our sign convention, D^D^Xq — D^Df^Xa = Riem^^a'^^/J ■ 
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of the metric. In an arbitrary local coordinate system, these quantities can be 
expressed as follows: 

[0.1. ZA) nicm^^^ — Ua^ ^ ^ Q !^ ^ Q M ^^lX^av' 

(3.1.2b) RiCp^ = Riem^^^" = daT^"^ - df^T^^ + T^^^T^^^ - T^^^TJ"^, 

(3.1.2c) i?''= 5"^Ric„^, 

(3.1.2d) r;^, ^.9"^5^5A. + d.g^x - dxg^,,). 

The twice-contracted Bianchi identities (see e.g. |Wal84] ) are 

(3.1.3) DaRic"'' - ^D'^i? = 0, (^i = 0, 1, 2, 3), 
which by p.l.ip leads to the following equations satisfied by T'^" : 

(3.1.4) DaT"f' = 0, (/X = 0,1,2,3). 

We will now briefly introduce the notion of a perfect fluid. For a more detailed 
introduction, including a very interesting account of the history of the subject, read- 
ers may consult Christodoulou's survey article |Chr07aj . The energy-momentum 
tensor for a perfect fluid is 



(3.1.5) T^'' = (p + 39)u"u'' + 

where p > is the proper energy density, p > is the pressure, and u is the 
four-velocity, a unit-length (i.e., UaU°' = —1) future-directed vectorfield on M.. The 
Euler equations, which are the laws of motion for a perfect relativistic fluid, are the 
four equations (|3.1.4p together with a conservation law p.l.6bp for the number of 
fluid elements. In an arbitrary local coordinate system, they can be expressed as 
follows: 



(3.1.6a) L)c,T"^ = 0, (^i = 0, 1, 2, 3), 

(3.1.6b) i:)„(nu") = 0, 

where n is the proper number density of the fluid elements. We also introduce the 
thermodynamic variable rj, the entropy per particle, which we will discuss below. 

Equations (|3.1.6ap - (j3.1.6b[) do not form a closed system, even in a fixed space- 
time [M.,gi^u)] there are too many fluid variables, and not enough equations. The 
standard way of remedying this difficulty is to appeal the laws of thermodynam- 
ics, which imply the following relationships between the fluid variables (see e.g. 
[GTZ99], IChrOTaj . |Chr07b! l: 

(1) p > is a function of n > and 77 > 0. 

(2) p > is defined by 

(3.1.7) P = n^ 

on 



- P. 

n 



where the notation |. indicates partial differentiation with • held constant. 
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(3.1.8) 



(3) A perfect fluid satisfies 
dp 



dn 



dp 



> with " = " ^=> Tj^O. 



As a consequence, we have that C-i the speed of soutkP^ in the fluid, is 
always real for 77 > : 

dp/dn\,, 



2 def- dp 



(3.1.9) C - 



dp/dn\r, 



> 0. 



(4) We also demand that the speed of sound is positive and less than or equal 
to the speed of light whenever n > and 77 > 0: 

(3.1.10) n>Oandr;>0 =^ 0<C<1- 

Postulates (1) - (3) are manifestations of the laws of thermodynamics and funda- 
mental thermodynamic assumptions, while Postulate (4) is connected to the notion 
of causality. More specifically, it ensures that at each x Ai, vectors that arc 
causal with respect to the sound con43 in T^Ad are necessarily causal with respect 
to the gravitational null con in TxA4. The physical interpretation of Postulate 
(4) is that the speed of sound is no greater than the speed of propagation of gravi- 
tational waves, i.e., that sound waves are causal as measured by g^,y. See |Spe09b| 
for a more detailed analysis of the geometry of the sound cone and the gravitational 
null cone. 

The assumptions p > 0, p > together imply that the energy-momentum tensor 
()3.1.5p satisfies both the weak energy condition {TapX'^X^ > holds whenever X 
is timelike and future-directed with respect to the gravitational null cone) and the 
strong energy condition {\Tfj^^ — l/2g°'^Tapg^i]X^X^ > holds whenever X is time- 
like and future-directed with respect to the gravitational null cone). Furthermore, 
if we assume that the equation of state is such that p = when p = 0, then p.l.9p 
and p.l.lOp together imply that p < p. This latter inequality implies that the dom- 
inant energy condition holds (— g^'^T^yX'^ is causal and future-directed whenever 
X is causal and future-directed with respect to the gravitational null cone). We 
remark that many important theorems in general relativity use the assumption that 
the matter model satisfies the dominant energy condition. As examples, we cite 
the Hawking- Penrose singularity theorem |Pen65j . |Haw67j and the positive mass 
theorem, which was first proved by Schoen-Yau [SY79|, |SY81j . and later by Witten 
|Wit81j . 

Under the remaining postulates. Postulate (1) is equivalent to making a choice 
of an equation of state, which is a function that expresses p in terms of 77 and p. An 
equation of state is not necessarily a fundamental law of nature, but can instead 
be an empirical relationship between the fiuid variables. In this article, we consider 



^■^In general, C, is not constant. However, for the equations of state we study in this article, 
is equal to the constant Cs ■ 

^^The sound cone is defined to be the subset of tangent vectors X £ TxM such that 

niaijX" = 0, where m^i^ g^i/ -|- (1 — C,^)uij,Ui, is the acoustical metric; see e.g. |Chr07b] . 
IRS09| . |Spe09b| for information concerning the role of the acoustical metric in the analysis of the 
Euler equations. 

^'^By gravitational null cone at x, we mean the subset of tangent vectors X S TxM such that 
go,pX°'XP = 0. 
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only the case of a barotropi^^ fluid under the equation of state p = c^p, 
< Cs < £ind according to (j3.1.9p , the constant Cg is the speed of sound. In 
this case, rj plays no role in our analysis of the Euler equations, and this quantity 
is absent from the remainder of the article. 

3.2. The initial value problem. In this section, wc discuss various aspects of 
the initial value problem for the Einstein equations, including the initial data and 
the notion of the maximal globally hyperbolic development of the data. We assume 
that wc are given an equation of state p = p{p) subject to the restrictions discussed 
in Section [31 We remark that the discussion in this section is very standard, and 
finds its in origin in the seminal works jCB52| by Choquet-Bruhat and |CBG69| by 
Choquet-Bruhat/Geroch. The discussion in this section concerns the initial value 
problem formulation "in the abstract," without reference to a choice of gauge for 
the Einstein equations. In Section [SJ we will provide a concrete formulation of the 
initial value problem for the Eulcr-Einstcin system that is suitable for our purposes 
at hand. Our formulation is based on a modification of the wave coordinate gauge 
used in [Rin08| , which is itself a modification of the gauge used by Choquet-Bruhat 
in |(]B52j . 

3.2.1. Summary of the Euler- Einstein system. We first summarize the results of the 
previous sections by stating that the Euler-Einstein system is the following system 
of equations: 

(3.2.1a) Ric - ^Rgf^i, + Ag^„ ^T^^, = 0, 1, 2, 3), 

(3.2.1b) 7^„r"^ = 0, (/x = 0,l,2,3), 

(3.2.1c) Dc,inu")=0, 

where T^^, — {p + p)u^Ui, + pg^u, and p — p{p) is given by an equation of state 
subject to the hypotheses discussed in Section IXTl 

We remark that by taking the trace of each side of p. 2. lap implies that i? — 4A = 
T, which implies that p. 2. lap is equivalent to 

(IMIal) Ric,,, - Ag^, - T,,, + ^Tg^, = 0. 

Since T = g°''^Tais = ip - p, it follows that 



(3.2.2) T^^ - 7;Tg^,v = {p + p)u^,u^ + - iP~P)9f^i^- 



]jTg^,„ = {p + p)Uf,u^ + i( 
Under the equation of state of interest to us, namely p = c^p, wc have that 



1 1 + 1 - 

(3.2.3) T^^ - -Tgf,^ = —^pu^u^ + ^^ Pg^u- 



16a 

barotropic fluid is one for which p is a function of p alone. 
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3.2.2. Summary and alternative formulation of the Euler- Einstein system under the 
equation of state p ^ c^p. Using the resuhs of Section [3.2.11 under the equation of 
state p = c^p, the Euler-Einstein system comprises the equations 

(3.2.4) Ric^. - A5^, - T^, + irg^, = 0, (//, = 0, 1, 2, 3), 

where T^^, ~ ^^^-^-^^u^u^+pg^y^ together with the equations of motion for a perfect 
fluid: 



(3.2.5a) DaT^^^O, (m = 0,1,2,3), 

(3.2.5b) Dc,{nu'^) = 0, 

and the normaUzation condition gapu'^u^ = — 1. 

By projecting in the direction of u and onto the 17— orthogonal complement of u, 
it can be checked that equations (|3.2.5aP are equivalent to the following system: 

(3.2.6a) u'^DaP + (1 + c^)pDaU°' = 0, 

(3.2.6b) (1 + cl)pu"Da,u^ + clWDa^p = 0, {j = 1, 2, 3), 

where 

(3.2.7a) gapu'^u^ = -1, 

(3.2.7b) n'""'^' u'^u'' + g"" . 

In the above expression, 11^'' is the projection onto the orthogonal complement 
of , which is sometimes referred to as the simultaneous space of u'^. Furthermore, 
for a barotropic equation of state, that is, one of the form p = p{p), the fundamental 
thermodynamic law p.l.7p . which reads p = — p in this case, can be used to 
show that equation p.2.5bp is an automatic consequence of p.2.5ap . Therefore, 
(|3.2.6ap - (j3.2.7bp is an equivalent formulation of the relativistic Euler equations 
under the equation of state p = c^p; we will work with this formulation for the 
remainder of the article. 

3.2.3. Initial data for the Euler-Einstein system. Initial data for the system p. 2. lap 
- (|3.2.1cp consist of a 3— dimensional manifold S together with the following fields 
on S : a Riemannian metric g^^ , a covariant two-tensor K_ji^ , a function p, and a 

vectorfield uP , {j,k = 1,2,3). 

It is well-known that one cannot consider arbitrary data for the Einstein equa- 
tions. The data are in fact subject to the following constraints, where D. is the Levi 
Civita connection corresponding to g : 

(3.2.8a) k - k^.ir'' + {r'K^bf = Spools, 

(3.2.8b) b^k^^ - tD.L., - To, Is- 

In the above expressions, indices are raised and lowered using g^^ and g. We remark 
that in the case of the equation of state p = c^p, and under the assumptions 
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discussed in the next paragraph, we have that Tools ~ ^^t^p(^o)^ ~Pj ^Qjlf. = 

^-^puQUj, where mq = — y^l + g^f^u^u'' ■ 

The constraints p.2.8ap - p.2.8bp are known as the Gauss and Codazzi equa- 
tions respectively. These equations relate the geometry of the ambient Lorentzian 
spacetime {Ai,g) (which has to be constructed) to the geometry inherited by an 
embedded Riemannian hypcrsurface (which will be (E, g) after construction). With- 
out providing the rather standard details (see e.g. |Chr08j . |Wal84| ) . we remark 
that equations p.2.8ap - (j3.2.8bp can be derived as consequences of the following 
assumptions: 

• S is a spacelikc submanifold of the spacetime manifold A4 

• g is the first fundamental form of E 

• ^ is the second fundamental form of E 

• The Eulcr-Einstcin equations are satisfied along E 

• We are using a coordinate system (a;° = t,x^ ,x'^ ,x^) on A4 such that 
Y. = {x e M \ t = 0}, and along E, goo = -1, goj = 0, gjk = g-i^, 
dtgjk = '^M^jkJ P = ^-iid = w-', with ga/3U°'u^ = —1. The ncxt-to- 
last condition is a consequence of the assumption that iu^,u^,u'^) is the 

orthogonal projection of the four-velocity {u'-', u^, u^, u^) onto E. 

We recall that under the above assumptions, g and K_ are defined by 

(3.2.9a) KUiX,Y)=gUiDj^X,Y), MX,YeT,t, 

(3.2.9b) gUX,Y) ^ gUX,Y) VX,reT,E, 

where TV is the future-directed normals to E at x, and D is the Levi-Civita con- 
nection corresponding to g. 

3.2.4. The definition of a solution to the Euler- Einstein system. In this section, 
we provide the definition of a solution to the Euler-Einstein system launched by a 
given initial data set. The discussion in this section and the next one follows the 
approach of [Rin08| , which was replicated in |RS09j . We begin with the following 
definition, which describes the maximal region in which a solution is determined 
by its values on a set S. 

Definition 3.2.1. Given any set S C A4, we define 'D{S), the Cauchy development 
of S, to be the union V{S) = V+{S) U V-{S), where V+{S) is the set of all points 
X £ M. such that every past inextendible causal curve through x intersects 5, and 
'D~ [S) is the set of all points p & M. such that every future- inextendibl^H causal 
curve through p intersects 5". 

Wc also define a Cauchy hypcrsurface. 

Definition 3.2.2. A Cauchy hypcrsurface in a Lorentzian manifold M is a. hypcr- 
surface E that is intersected exactly once by every inextendible timelike curve in 
M. 

^'''Under the above assumptions, it follows that at every point a; S E, TV^ = 5^. 

^^A curve 7 : [so, Smax) — > M is said to be future-inextendible if there does not exist an 
immersed future-directed curve -y : I M with [soi^jnua:) C /, [so,Smax) ^ I, and 7|[so,s,„„^) = 
7. Past inextendibility is defined in an analogous manner. 
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It is well-known that if S C is a Cauchy hypersurface, then = (see 

e.g. lO'NSSQ . 

We now provide the definition a solution launched by an initial data set. 
Definition 3.2.3. Given sufficiently smooth initial data (S, g , (j, k = 

1,2,3) as described in Section [3.2.31 a (classical) solution to the system p. 2. lap - 
p.2.1cp is a 4— dimensional manifold A^, a Lorentzian metric g^^, a function p, a 
vcctorficld w^, (/i, = 0, 1, 2, 3), and an embedding S ^ subject to the following 
conditions: 

• (7 is a tensorfield, p is a function, and u is a tensorfield 

• Equations l|3.2.1ap - p.2.1cp are satisfied by the components of 5, p, and u 

• S is a spacclike Cauchy hypersurface in [M , g) 

• g IS the first fundamental form of S 

• ^ is the second fundamental form of E 

• 7r(u) = u where tt denotes g— orthogonal projection onto E 

The array (A^, g^^jp, w^) is called a globally hyperbolic development of the initial 
data. 

3.2.5. The maximal globally hyperbolic development. In this section, we state a fun- 
damental abstract existence result of Choquet-Bruhat and Geroch |CBG69| , which 
states that for initial data of sufficient regularity, there is a unique "largest" space- 
time determined uniquely by it. The following definition captures the notion of this 
"largest" spacetime. 

Definition 3.2.4. Given sufficiently smooth initial data for the Eulcr-Einstein 
system (|3.2.1ap - p.2.1cp that satisfy the constraints (|3.2.8ap - (j3.2.8bp . a maximal 
globally hyperbolic development of the data is a globally hyperbolic development 
{A4,g,p, u) together with an embedding l : T, ^ Ai with the following property: if 
(A^', g',p', u') is any other globally hyperbolic development of the same data with 
embedding t' : S — > A^', then there is a map V' : A^' A^ that is a diffeomorphism 
onto its image such that i]j*g = g','ijj*p = p','ip*u = u' and ijj o l' = l. Here, ^* 
denotes the puUback by i/'. 

Before we can state the theorem, we also need the following definition, which 
captures the notion of having two different representations of the same spacetime. 

Definition 3.2.5. The developments {Ai, g,p,u) and {A4' , g' ,p' ,u') are said to be 
isometrically isomorphic if the map if) from the previous definition is a diffeomor- 
phism from J\A to M.' . 

We now state the aforementioned theorem. 

Theorem 3.1 (Existence of an MGHD). jCBG69j Given sufficiently smootl^ 
initial data for the Euler- Einstein system p. 2. lap - p.2.1cp . there exists a max- 
imal globally hyperbolic development of the data which is unique up to isometric 
isomorphism. 

The remainder of this article concerns the "future" properties of the maximal 
globally hyperbolic developments of sufficiently smooth data near those correspond- 
ing to the FLRW background solutions introduced in Section 21 

^'^The article ICBG691 only discusses the case of smooth data. However, as discussed in 
ICGPlOl Section 6] , the regularity assumptions of Theorem 15.11 are sufficient for the conclusions 
of Theorem 13. ll to be valid. 
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4. FLRW Background Solutions 

Our main results address the future-stability of a well-known class of FLKVV0 
background solutions ([0, cxd) xT'^, u^), (/i, i/ = 0, 1, 2, 3), to the system p.2.4|) 

+ p.2.6a|) - (j3.2.6b|) . These background solutions physically represent the evolution 
of an initially uniform quiet fluid in a spacetime that is undergoing exponentially 
accelerated expansion. To find these solutions, we follow a standard procedure that 
is outlined e.g. in jWal84[ Chapter 5] which, under appropriate ansatzes, reduces 
the Euler-Einstein equations to ODEs. The discussion in this section was essentially 
provided in jRS09[ Sect ion 4], but we repeat it here for convenience. To proceed, 
we first make the ansatz that the background metric ]j = g{t) is of the form 

3 

(4.0.10) g = -dt'^ + a^{t)^{dx'f, 

i=l 

from which it follows that the only corresponding non-zero Christoffel symbols are 

(4.0.11) fj\=f,'j=adS,k, r/o = To^-^. = (j,fc = l,2,3), 

where 



dcf a 
w = -, 
a 

Using definitions (|3.1.2bp and (|3.1.2cp . together with (|4.0.1ip . we 

(4.0.13a) Rkoo- J^5oo = 3f-)^ 

(4.0.13b) mcoj - ^Rgoj = 0, (j = 1, 2, 3), 

(4.0.13c) mc,k-^Rgjk = -{2ad + d^)Sjk, (j,fc = 1,2,3). 

We then assume that p = p{t), p = pit), and u^ = (1, 0, 0, 0). We also assume that 
the equation of state p = c^p holds, and for simplicity, we assume for the remainder 
of the article that 



(4.0.12) 
and a = -±a. 

at 

compute that 



(4.0.14) a(0) = 1. 

Inserting these ansatzes into the Einstein equations p.2.4|) . we deduce (as in [Wal84| ) 
the following equations, which are known as the Friedmann equations in the cos- 
mology literature: 



(4.0.15a) pa3(i+'=^) = p. 



(4.0.15b) 



3 3 V 3 3a3( 



l+c- 



2\ I 



^^Technically, the term "FLRW" is usually reserved for a class of solutions that have spatial 
slices diffeomorphic to K.^, or hyperbolic space (see |Wal84| ) . 
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where the positive constant p denotes the initial (uniform) energy density. We 
also denote the initial pressure by p == c^/j. Observe that the rapid expansion of 
the background spacetime can be easily deduced from the ODE ()4.0.15bj) . which 
suggests that the asymptotic behavior a{t) ^ e^*,H =^ A more detailed 

analysis of a(t) is given in Lemma 14.1.11 

For aesthetic reasons, we also introduce the quantities 



(4.0.16) n{t)=lnait), 
(4.0.17) >^ = 3c2, 

which implies that 

(4.0.18) a(0 = e"(*), 

(4.0.19) u{t) = j^n{t). 

For future use, we note the following simple consequences of the above discussion; 
we leave it to the reader to supply the details: 

(4.0.20a) 3a;2 - A = ^p, 

(4.0.20b) uj = -^f^P, 

(4.0.20c) 3^ + 3cj2_^^_l±|£|~ 

2ci 

4.1. Analysis of Friedmann's equation. The following lemma summarizes the 
asymptotic behavior of solutions to the ODE (j4.0.15b[) . 

Lemma 4.1.1. jRS09[ Lemma 4.2.1] Let p,(^ > be constants, and let a{t) be the 
solution to the following ODE: 



(4.1.1) l-Vt + a!^' «('^) = ^- 

Then with H = -^/A/S, the solution a{t) is given by 



2/^ 



(4.1.2) yXTT,„,.(£|?)} 

and for all integers N > 0, there exists a constant Cn > such that for all t > 0, 
with A=i^^ (\/ 3^ ^ ^ ^-^ } ' have that 

(4.1.3a) (l/2)2Ae^* < a{t) < Ae"\ 

(4.1.3b) e-"'^a{t)- AH^ < Cnc-'^K 
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Furthermore, for all integers N > 0, there exists a constant Cjv > such that 
for all t > 0, with 



(4.1.4) 

we have that 



def a 



(4.1.5a) 
(4.1.5b) 



H<uj{t)<\lm + !^, 



< Cnc 



-<;Ht 



Remark 4.1.1. Because of equation (|4.0.15b[) . we will assume for the remainder 
of the article that ^ = 3(1 + c^). 



5. The Modified Euler-Einstein System 

In this section, we recall the wave coordinate system introduced in jRS09| , which 
was based on the framework used in [RinOSj . We then use algebraic identities 
that are valid in wave coordinates to construct a modified version of the Euler- 
Einstein equations, which is a system of quasilinear wave equations coupled to 
the first-order Euler equations, and which contains energy- dissipative terms. We 
then construct data for the modified system from given data for the unmodified 
system in a manner consistent with our wave coordinate system. Next, to facilitate 
our analysis in later sections, we algebraically decompose the modified system into 
principal terms and error terms. We then discuss local existence and a continuation 
principle for the modified system, and we sketch a standard proof of the fact that the 
modified system is equivalent to the unmodified system if the Einstein constraint 
equations and the wave coordinate condition are both satisfied along the initial 
Cauchy hypersurface S. Finally, for convenience, we introduce some PDE matrix- 
vector notation for the Euler equations. 

5.1. Wave coordinates. To hyperbolize the Einstein equations, we use a coor- 
dinate system in which the contracted ChristofFcl symbols =^ 9°'^^a p 

spacetime metric g are equal to the contracted ChristofFcl symbols F'' g"^^^ p 
of the FLRW background metric g. This well-known condition is known as a wave 
coordinate condition since F^ = F^ if and only if the coordinate functiontl^ a;^ are 
solutions to the wave equation g'^^DaDpx'^ f = 0. Using (|4.0.10p and (|4.0.1ip . 
we compute that in wave coordinates, we have 

(5.1.1) F^ = F^ = 30.5^, F,, = g^„F" = 3ugo^, 

where uj{t) is defined in (|4.0.12|) . 



The are scalar- valued functions, despite the fact that they have indices. 
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We now introduce the tensorfielda i 

(5.1.2a) = 3uj6^, P„ = r„ = 8^30^, 

(5.1.2b) = - ^^ = - Tp. 

The idea beliind wave coordinates is that when = 0, whenever it is expedient, 
we may replace F'* with 3uj5q (and vice-versa) without altering the content of the 
Einstein equations. The existence of such a coordinate system is nontrivial, and 
it was only in 1952 that Choquet-Bruhat jCB52| first showed that they exist in 
general (see Proposition I5.6.1"|) . With this idea in mind, we define (as in |Rin08[ 
Equation (47)]) the modified Ricci tensor Ric^j/ by 

(5.1.3) 



where 



(5.1.4) Ug"^ g'^^d^dp 

is the reduced wave operator corresponding to the metric g. 

We now replace the Ric^ii/ with Ric^,y in p.2.4p . expand the covariant differen- 
tiation in p.2.6ap - (j3.2.6bp . and add additional inhomogencous terms /^j^ to the 
left-hand side of p.2.4p . thereby arriving at the modified Euler- Einstein system: 

(5.1.5a) Ric^ii/ - Ag^^ - T,,,^ + ^Tg^,, + /^^ = 0, (//, = 0, 1, 2, 3), 
(5.1.5b) u'^Do,p+{l + cl)pDau" ^0, 

(5.1.5c) u-^D^u' + ' n^'^D^p = 0, (j = 1,2,3). 

(1 + ci)p 

Here, the additional terms arc defined to be 

(5.1.6a) /oo = -2wQ° = 2w(r" - 3a;), 

(5.1.6b) Iqj = Ijo 2ujQj = 2uj{3ugoj - T^), 

(5.1.6c) I^k^I^^'^'o. 

We have several important remarks to make concerning the modified system 
(|5.1.5ap - (|5.1.5cp . First, because the principal term on the left-hand side of (|5.1.5ap 
is — iQgg^;^, the modified equations (j5.1.5ap are quasilinear wave equations, and are 
therefore of hyperbolic character. Since the Euler equations (|5.1.5bp - (j5.1.5cp are 
also hyperbolic in a fixed spacetime, (see e.g. }Chr07b| . Spe09b| ), it follows that 



^■^Technically, and do not have the coordinate transformation properties of a vcctor- 
field/covectorfield. Nonetheless, we will treat them as vectorfields/covectorfields when we compute 
their covariant derivatives. On the other hand, Q^" and do have the transformation properties 
of a vcctorfield/covectorfield. 
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the modified Euler-Einstein system is a hyperbolic system of mixed order. Second, 
the gauge terms /^j/ have been added to the system in order to produce an energy 
dissipation effect that is analogous to the effect created by the 3{dtv)^ term on the 
right-hand side of the model equation (|1.2.ip . These dissipation-inducing terms 
play a key role in the global existence theorem of Section [TTJ Finally, in Section 
15.51 we will elaborate upon the following fact: if the initial data satisfy the Gauss 
and Codazzi constraints p.2.8ap - ()3.2.8bp . and if the wave coordinate condition 
Q^\t=Q = is satisfied, then Q^, /^j^ = 0, and Ric,,^ = Ric^j^; i.e., under these 
conditions, the solution to (|5.1.5al) - (|5.1.5cP is also a solution to the Eulcr-Einstcin 
system (|3.2.1ap - (|3.2.1c|) . 

5.2. Summary of the modified system for the equation of state p = c^p. 

For convenience, we summarize the results of the previous section by listing the 
modified Euler-Einstein system (where j = 1, 2, 3): 

(5.2.1a) Ricoo + 2a;r° - 6uj^ 

(5.2.1b) Ricoj - 2w(rj - Sujgoj) 

(5.2.1c) mcjk 
(5.2. Id) 
(5.2. le) 

where gapu'^u^ = -1 and n^"" = u^'u'' + g^". 

5.3. Construction of initial data for the modified system. In this section, 
we assume that we are given initial data (S, p_^.^, i;^^^,,^, u-'), (j, fc = 1,2,3), for 
the Euler-Einstein equations p. 2. lap - (|3.2.1cp as described in Section 13.2.31 In 
particular, we assume that they satisfy the constraints p.2.8ap - (|3.2.8bp . We will 
use this data to construct initial data for the modified equations that lead to a 
solution (Al, u^), (//, = 0,1,2,3), of both the modified system and the 
Einstein equations; recall that a solution solves both systems <==^ = 0, where 

is defined in (j5.1.2bp . We remark that one may consider arbitrary data for 
the modified equations (|5.2.1ap - (|5.2.1cp . but without further assumptions, the 
resulting solution is not necessarily a solution to the Einstein equations p. 2. lap - 
(ISTO)) . 

To supply complete data for the modified equations, we can specify along S = 
{t = 0} the full spacetime metric components (7^jy|t=o, their future-directed normal 
derivatives dtgfj,iy\t=o, the pressure p, and the 17— orthogonal projection of the four- 
velocity onto E. To satisfy the requirements 

• ±^{t = 0} 

• g is the first fundamental form of S 

• _K is the second fundamental form of S 

• dt is future-directed and normal to S 

• P\±=P 



- Affoo -p{^—^{uof + -^^500 j = 0, 

/ 1 + 1 - \ 

- Agjk - p{^—^UjUk + 2c^' 9]kj = 0, 

w"i:'„p+(l + c^)pD„u" =0, 
u'^D^u^ + f\ WD^p = 0, 
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• u is the 5— orthogonal projection of (the unit-normahzed) four-velocity u 
onto S, 

we set 

(5-3. la) 5oo|t=o = -1, 5oj |t=o = 0, g^^lt^o ^ g^,., 

(5.3.1b) p\t=o=P, Uj\t=o = Uj, Mo|t=o = --y/l +°g^Jku', 

(5.3.1c) {dtg,k)\t=o = 24.,. 

To satisfy the initial wave coordinate condition Qfj\t=o = 0, we first compute 
that 

(5.3.2a) ToUo = -\{dtgoo)\t=o ~ t^ab. 

(5.3.2b) r. li.o = -{dt9oi)\t=o + \ti'^dag,^ - 9 flab)- 

Using (|5.3.2ap and (j5.3.2b[) . the condition Qfi\t=o = is easily seen to be equivalent 
to the following relations, where w is defined in (|4.0.12p : 

(5.3.3a) 

{dt9m)\t=o = 2(-3^|t=o5ook=o -f'^aJ = 2(3c.(0) - f^,,), 

-1 

(5.3.3b) 



We remark that in the above expressions, g^'^ denotes a component of the inverse 
of g. This completes our specification of the data for the modified equations. 

5.4. Decomposition of the modified system in wave coordinates. In order 
to reveal the dissipative structure discussed in Section [TT^ we decompose the mod- 
ified equations (|5.2.1ap - (|5.2.1eP into principal terms and error terms, which we 
denote by variations of the symbol A. A central component of our global existence 
argument is the derivation of suitable bounds for the error terms; the estimates 
of Section |9] will justify the claim that the A terms are in fact error terms. We 
begin by recalling the previously mentioned rescaling hjk of the spatial indices of 
the metric: 

(5.4.1) h,,'='e''%,. 

We will also make use of the following rescaling of the pressure: 

(5.4.2) P e3(i+^')^V- 

The above rescaled quantities will be order 1 in our global existence theorem, which 
makes them more convenient to work with; i.e., we have rescaled purely for conve- 
nience. We remark that for the FLRW background solution, P = p and hjk = Sjk- 
The decomposition is carried out in the next proposition. 
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Proposition 5.4.1 (Decomposition of the Modified Equations). The equa- 
tions (|5.2.1a[) - (j5.2.1cp in the unknowns {g^j,^^ P, u-'), {fi, v ^ 0,1, 2, 3), {j = I, 2, 3), 
can be written as 



(5.4.3a) Dgigoo + 1) = 5Hdtgoo + dH^goo + 1) + Aqo, 

(5.4.3b) Dggoj = SHdtgoj + ^H^go, - 2Hg-'Tajb + ^Oj, 

(5.4.3c) Dghjk = iHdthjk + Ajk , 



(5.4.3d) 

u"a„(P - p) + (1 + c2) f — ^PUadtu"- + (1 + C^)PdaU'' 



A, 



(5.4.3e) u"da,u' + '-—W"d^iP - p) ^ (x - 2)uju^ + A', 

(1 + cf)P 

where 



(5.4.4) = + Jl + ('-^Y - (500 + ly 

500 V ^ 500 ^ .goo ^ .goo ^ 

(5.4.5) U^"" ^ u^'u•' + g^"" , 



(5.4.6) ^=Vf' --3c^' 

the error terms A^^, A, A^ can be expressed as 

(5.4.7a) 

Iaoo = A^,oo + Ac.oo - ^^(300 + l)e-3(i+c^)np- _ M±I^_3(i+c^)0(p _ -^^ 

- i±^(uo + 1)(«„ - l)e-^(i+^=)"P - ^(500 + l)e-3(i+^^)"P 
+ |(l. - H)dtgoo + Huj^ - H^)igoo + 1), 

(5.4.7b) 

- i^e-'^'+^'^''Pffo, + - H)dtgoj + (^' - i?').go, - (^ - i^),?"^,,;., 

(5.4.7c) 

\a,u = e-2^^A^,,., + i±fipe-3(i+^^)^^(500 + l)h,^ - 2c.g°'^a,/.,, 

- ^e-3(i+^=)^^(P -p-)/,,., - i±ie-2a^-3(i+cj)np^^.^^ ^ _ ^^^^^ 
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(5.4.7d) A = -(1 + c2)FA„V° - (1 + cl)PA^\u'^ 



2uq 



-{(at5oo)(«")' + 2(atgoa>%" + (at.ga6)w"w^}, 



(5.4.7e) A^' = (>f - 2)w(u° - l)?/^' - pu'^u^ + >fa;.g°^ 

</ie A^.^^ are defined in (j5.4.17ap - (j5.4.17cp . Ac.oOj '^COj fJ^'e defined in (j5.4.19ap 
- (|5.4.19bp . and the are defined in (|5.4.20ap - (j5.4.20f|) . 

Furthermore, vP is a solution to the following equation: 



(5.4.8) 

where 



(l + c2)P 



n""a„(p-p) = A'J 



(5.4.9) A" = >^cj{(u" - l)(u" + 1) + (ffO" + 1)} - iogabu'^u'' - A„%?/"u'5. 

Remark 5.4.1. Equation (|5.4.4p enforces the normalization condition gapu^u^ = 
— 1 and the future-directed condition u° > 0. 

Proof. The proof is a series of tedious computations based Lemma 15.4.21 - Lemma 
15.4.51 We provide the proofs of (|5.4.3cp - (|5.4.3ep and leave the remaining details 
to the reader. To obtain (|5.4.3cp . we first use equation (|5.2.1cp . Lemma [5.4.21 and 
Lemma [5.4.31 to obtain the following equation for hjk = e^^^^gjk ■ 



(5.4.10) 

Dghjk = Siodthjk + 2[3a;2 + dtcu - A]hjk - ^Log^^-dah^k 

+ 2e-^/\j,,,k - 2(9ta;)(ff"" + l)h,k - ^ii±£lle-2^^e-3(i+^^)^^P«,«, 

Now using (|4.0.20ap and (|4.0.20bp . it follows that 2[Zuj'^+dtOJ-K\hjk = ^-^pe-'^^^+^^^^hjk- 

Substituting into (|5.4.10l) . and using dtw = -ii±§llpe-3(i+=')" (i.e., (I4.0.20bp ). it 
follows that 



(5.4.11) Ugh.k = iujdth.k 2^e~3(i+^''^'(P -p)/i,fe _ ^^gOo.Q^^^^ 

+ 2e-'''AA;,k + ii±^p-e-3(i+^=)^^(g°" + l)h,k 

_ 2(l + c|)^_2^^_3(i+,.)f,^^^^^^ 

Equation (|5.4.3cp now easily follows from (|5.4.11l) . We remark that the proofs of 
(|5.4.3ap and (|5.4.3bp require the use of Lemma 15.4.41 

To obtain (j5.4.3dp . we first expand the covariant differentiation in (j5.2.1dp to 
deduce the following equation: 
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(5.4.12) u"d^p+{l + cl)pd^u'' = -{l + ci)pr^^pu^. 

Lemma 15.4.51 impliGS that — iuJvP + A^^u^ , while the normalization con- 

dition gaj3u"u^ = —1 imphes that dtu^ = —-^{uadtu"' + ^{dtgai3)u"u^} . Equation 
()5.4.3dp now follows from multiplying both sides of ()5.4.12p by e^^^^^^-*^ and using 
the fact that e^^^+'^i^^ dtp = dtP - 3(1 + c'i)ujP. 

Similarly, to obtain (j5.4.3e|) . we first expand the covariant differentiation in 
(j5.2.1cp to deduce the following equation: 

(5.4.13) u'^d^u^ + j^r^^^'"d^p = -ri/3u"?/^. 

Lemma [5X5] implies that Lj pu"u^ = 2oju°u^ + A J pU°'u'^ . Equation (|5.4.3ep now 
follows from (|5.4.13p and the fact that (■i_|_e2)p n^"^aP = ji^:^yp^''°'do,{P — p) — 

□ 



We now state the following four lemmas, which are needed for the proof of 
Proposition 15.4.11 

Lemma 5.4.2. |Rin08| Lemma 4] The modified Ricci tensor from ()5.1.3p can be 

decomposed as follows: 

(5.4.14) 

1 _ 3 3 

Ric^L' = --^Oggfji, + -{go^d^uj + goi,d^uj) + -oodigf^u + Afj,^, (//, = 0, 1, 2, 3), 
where 
(5.4.15) 

= g^^g"^ [{daguK){dpg^x) - r,,,r^^A] , (m, = o, i, 2, 3). 

□ 

Lemma 5.4.3. |Rin08[ Lemma 5] The term A^i, (ji^v = 0, 1,2,3J defined in 
(|5.4.15p can be decomposed into principal terms and error terms A^.^^ as follows: 



(5.4.16a) 

Aqo = 3a;2 - ug'^'^dtgab + 2cj.g"''9a.go6 + Aa,oo, 
(5.4.16b) 

= 2wg°°at5oj - aw^gOO^oj - w5°°^j5oo + ws^'r,,!, + A^.o^, [j = 1, 2, 3), 
(5.4.16c) 

A,k = 2wg™9t.g,fe - '2^''g°''gjk + A^,,fe, (j, fc = 1, 2, 3), 
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where 



(5.4.17a) 

Aa,oo = ig°y{{dt9oof - (Fooo)'} + .g°°.g°°{2(at5oo)(atffoa + 9a5oo) - 4rooorooa} 
+ 5°°5"'{(5t<?oa)(atgo6) + {dagooKdbgoo) - 2rooaroo6} 
+ g°''g°''[2{dtgoo){dagob) + 2{dtgob){dagoo) - 2roooraob - 2rno6rooa} 
+ g"'''g°'[2idtgoa)idigob) + 2{dbgm){dagoi) - ^TmaTmb] 

+ 5"''5'™(5a5o/)(45orn) + g^'Otgai - 2c^<S? ){dbgom + dragob) 

2 ^ ' 

- ^5°''5''"(5a3o; + digoa){dbgom + 9„i5oh) 

-7( g''%gai-2u:S^ )( g'^^dtgbm. ~2iodi ), 
4 „ ' V ^ ' 



(5.4.17b) 

Aa,oj = ig°°y{{dtgoo)idtgoj) - Tooorojo} 

+ ff°°5°''{(5t5oo)(9tgQj + dagoj) + idtgoj)idtgoa + dagoo) 

~ 2roooroja — 2rojorooa| 

+ g°°{g^%gbj-2uj5])(^dtgoa - ^^affoo) + lg°°g''\dagoo){dbgoj + d^gob) 
+ g"''g°''[{dtgoo){dagbj) + {dtgob)idagoj) + {dagoo)idtgbj) + {dagob){dtgoj) 

+ g"'' g°'' {{dtgoa){digbj) + idigoa)idtgbj) + {dbgm){dagij) + (dbgaOidagoj) - 2rooar/j6| 

- g"'^g°'^[{digoa + dagoi)^ojb - ^{dtgia){dbgQj - djgob)^ 

+ ujg°-{dtgaj - 2Logaj) + iff"'( g'^'dtgia - 2ujSf )dtgbj 
V ' 2 V ^ ' 

e^"dthaj e^-^'g-"'dth,^-2u:gOi'go, 

+ <?'^''<?""{(Sa50i)(565mj) - lidagol + digoa)Tbjrn} + ( /"^tgig - 2^.^ )r,,^ , 

e^^ g'"^ dthia— g°"^ goa 
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(5.4.17c) 

+ 9"°g°°'[{dtgoj)idt9ak + dagok) + {dtgok){dtga] + dagoj) 

^ 2rojorofca — 2rofeoroj"a| 

+ 9'^°g°'''\{dago3){dbgok) - ^(^aSoj - djgna){dbgok - 5fe.9oh)} 

g'^'dtgao-'^^S) ){digok - dkgob) + { g'^'dtg^u - 2uj51 )idago, - 
^ ^ ^ V ' ^ « ' 

+ c.ff°°( gfcfcg°^^-(5g )at.gaj- + ^g"°( g°^9tga, - 2c.,5^^ )idtgbk_^^u^) 

-gokg"" e^^g''''dthaj-2Ljg'>'>goj e'^^^dthbk 

+ ff°°S°''{(5tgoj)(9a.gbfc) + {dtgbj){dagok) + idagoj){dtgbk) + idagb3){dtgok) 
+ g''''g°'^[idtga])idigbk) + {diga]){dtgbk) + {dbga]){dagik) + {dbgi]){dagok) 

— 2rojar(fc6 — 2r;jarofc6| 
+ g"^^ g""'^ ^{dagij)(dbgmk) - TajlTtkm^ 

□ 

Lemma 5.4.4. |Rin081 Lemma 6] The sums Aqq + Iqq and Aqj + loj, {j ~ 1, 2, 3), 
can be decomposed into principal terms and error terms as follows, where Ioo,Ioj 
are defined in (|5.1.6ap - (|5.1.6bp ; Aoo^Aoj are defined in (j5.4.15p .' and Aa,oo, '^A,oj 
are defined in (|5.4.17a|) - (|5.4.17b|) ; 

(5.4.18a) ^00 + 2cjr" - Gcj^ = Lodtgoo + 3uj'^{goo + 1) + Sw^goo + Aa,oo + Ac,oo, 
(5.4.18b) 

Aoj + 2uj{3uigoj - Tj) = 4a;^goj - ujg^-^Tajb + Aa.oj + Ac,oj, 

where 
(5.4.19a) 

Ac,oo = -6(5oo)"'^'{(5oo + l)'-5%oa}-^(g™ + l)( g'^'dtgab ^ 6u ) 

e^^'g'-i-dthab-2ujgO'^goa 

+ 2c^(.g™ + l)g'''dagob + ^(.9™ + 1)(<?™ - l)5tffoo + 2a;g™g°'' (Loao + 2Tooa) 
+ 4c.50"g°^roa, + 2c^ff"''5°'r,z,, 
(5.4.19b) 

Aaoj - 2^2(^00 ^ _ 2cjg"°|( 9tg„j- -2c^ffajO + dagaj - d^goa}- 
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□ 

Lemma 5.4.5. The C'hristojfel symbols can be decomposed into principal 

terms and error terms A^"^ as follows: 

(5.4.20a) = Ao%, 

(5.4.20b) r/„ r„o^. = A/o = Ao", , 

(5.4.20c) r,\ = A,\, 

(5-4.20d) = r,^; = uSi + A,\ = u5i + „ 

(5.4.20e) r/, = + A/„ 

(5.4.20f) r,^, = A>„ 



where 



.g°"at.goo + 2g°''dtgoa - .g°°9a.goo, 

g'^dtgm + 2g'''dtgoa-g"'dagm, 

g'°dkgoo + .9^° 9,50a - g'^'dagok + ( 5'"9tffafc - 2u;Si ), 

^ ^ ^ 

9°°{djgok + dkgoj) + g°"-{djgak + dugaj - dagjk) 
+ {dtgjk ' 2u;g,k) - 2(5°° + l)ujg,k - (5"" + iKdtgjk - ^Log^k), 

" V ' " V ' 

(5.4.21f) 

2A,';. = + d,g,a ~ ^ag^J). 

Proof. The proof is again a series of tedious computations that follow from the 
definition T^^ = ^"^{d^gxty + d^gf^x - 9a5m"^)- ^ 

5.5. Classical local existence and the continuation principle. In this section, 
we discuss classical local existence results and continuation criteria for the modified 
system (|5.4.3ap - (|5.4.3ep . The theorems in this section are stated without proof; 
we provide references for the rather standard techniques that can be used to prove 
them. 

Theorem 5.1 (Local Existence for the Modified System). Let N > 3 be an 

integer. Let g^^ = g^,„\t=o, 2K^^ = {dtgfiiy)\t=o, {p-^v = 0,1,2,3), P = P|t=o = 
p|t=Oj — u^t=o, {j — 1,2,3), gapu"u^\t=o = —1, be initial data (not necessarily 



(5.4.21a) 

2A ° - 
(5.4.21b) 

2A ° - 

(5.4.21c) 

2A ^ — 

(5.4.21d) 

2A ^ — 

(5.4.21e) 

2AA = 
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satisfying the Einstein constraints) on the manifold for the modified system 
(|5.2.1ap - (|5.2.1cp satisfying (for j,k ^ l,2,'S) 

(5.5.1a) dgjk 6 H^" , 500 + 1 e H''+\ goj G H^+\ 

(5.5.1b) K,k - L^m^k e H^, Koo e i?^, Koj € if^, 

(5.5.1c) P ~peH^ , , 

where p > is a constant. Assume that mij-^ja p > 0. Assume further that there is 
a constant C > such that 

(5.5.2) cSabX^x'' < ijabX^x^ < c~'^SabX''x'', y{x\x^,x^) e m.^, 

and such that goo < 0, so that by Lemma \7.0.1\ the 4x4 matrix is Lorentzian. 
Then these data launch a unique classical solution {gfj^^,P,u^) to the modified sys- 
tem existing on a slab [T-,T+] X T^, with T_ < < r+, such that 

(5.5.3) 

.9p. € C2([r_,r+] X t3), p e c\[t^,t+] x t^), e c\[t^,t+] x t^), 

such that goo < 0, and such that the eigenvalues of the 3x3 matrix gjk are uniformly 
bounded below from and from above. 

The solution has the following regularity properties: 

(5.5.4a) 

dg,k e C°([T_, T+], i?^), goo + le C"([T_, i/^+i), go, £ C"([T_, r+], i/"+i), 

(5.5.4b) 

dtgjk - Mt)93k e C°([r_,r+],i7^), dtgoo G C"([T_,r+],iJ^), 5tgo, e C"([T_,T+],ff^), 
(5.5.4c) 

P-pe C°([r_,T+],if^), 6 C"([T_,T+],if^), 7/0-1 eC°([T_,T+],ff^). 
Furthermore, g^i, is a smooth Lorentzian metric on {T-,T^) x T"^, and the sets 

{t} X are Cauchy hypersurfaces in the Lorentzian manifold {A4 = {T^^T^) x 
T\g^,) forte {T^,T+). 

In addition, there exists an open neighborhood O 0/ (g^i,, ilT^j^, P, ■«■') such that 
all data belonging to O launch solutions that also exist on the slab [T_ , T+j x 
and that have the same regularity properties as (.g^,/, P, u^). Furthermore, on O, the 
map from the initial data to the solution is continuous^ 

Finally, if, as described in Section \5.SX the data for the modified system are con- 
structed from data for the Einstein- Euler system satisfying the constraints p.2.8ap 
- ()3.2.8bp on an open subset 5 G T'^, and if the wave coordinate condition Q^\s = 
holds, then {gfj^^,p,u^) is also a solution to the unmodified equations (j3.2.4p - 
(|3.2.5bp onV{S), the Cauchy development of S. 

Remark 5.5.1. The hypotheses in Theorem [511] have been stated in a manner that 
allows us to apply to it initial data near that of the background solution of Section 
m Furthermore, we remark that the assumptions and conclusions concerning the 



^■^By continuous, we mean continuous relative to the norms on the data and the norms on the 
solution that are stated in the hypotheses and above conclusions of the theorem. 
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metric components gjk would appear more natural if expressed in terms of the 

variables hjk e~'^^gjk', these rescaled quantities are the ones that we use in our 
global existence proof. 

Proof. The existence aspect of Theorem [O] can be proved using standard methods 
that follow from energy estimates in the spirit of those proved below in Sections 
01131 and^ See e .g. |H6r97l Ch. VI], |Maj84[ Ch. 2], [SS981 Ch. 5], |Sog08[ 
Ch. 1], |Spe09b| , and |Tay97[ Ch. 16] for details on how to prove local existence 
as a consequence of the availability of these kinds of energy estimates. Also see 
|Rin08[ Proposition 1]. The fact that u'"), where p = e^^(^+''-')^^P, is also 

a solution to the modified equations if the constraints and the wave coordinate 
condition Q^|s ~ arc satisfied is discussed more fully in Section [5.61 

Remark 5.5.2. Our use of energy currents in Sections 16.31 and [TUl to derive energy 
estimates for the relativistic Euler equations may be unfamiliar to some readers. 
However, once one has such estimates, local existence can be proved using the 
standard arguments mentioned above. 

□ 

In our proof of Theorem I ll.li we will use the following continuation principle, 
which provides standard criteria that are sufficient to ensure that a solution to the 
modified equations exists globally in time. 

Theorem 5.2 (Continuation Principle). Assume the hypotheses of Theorem 
\5.1[ Let Tmax be the supremum over all times such that the solution {g^^^^pjU^) 
exists on the interval [0,T_|-) and has the properties stated in the conclusions of 
Theorem \5.1i Then ifTmax < oo, one of the following four possibilities must occur: 

(1) There is a sequence {tn, a;„) S [0, T^ax) xT'^ such that lim„^oo <7oo(in, a^n) = 
0. 

(2) There is a sequence {tn,x„) S [0,Tmax) x TT'^ such that the smallest eigen- 
value of gjk{tn,Xn) converges to as n co. 

(3) There is a sequence {tn, Xn) € [0, T„iax) x T'^ such that lim„_^oo P{tm Xn) = 
0. 

(4) lim^^y-^^ supo<^<4 I E|i,^=o {\\9i^'^{'r, ^Wcf + \\dtg,,u{T, ■)\\cl) 

+ \\P{t, ■)\\ci + \\dtP{r, OIU^+E ■=! {\\uKr, ■)\\cl + \\dtu^ {r, | = oo. 

Similar results hold for an interval of the form {Tmin , 0] . 

Proof See e.g. |Hor97l Ch. VI], |Sog08[ Ch. 1], |Spe09a| for the ideas behind a 
proof. If either (1) or (2) occurs, then the hyperbolicity of the operator Og can 
break down. Condition (3) is connected to the fact that the Euler equations can 
degenerate when the pressure vanishes. □ 

5.6. Preservation of the wave coordinate condition. In Section 15.31 from 
given initial data for the Einstein equations, we constructed initial data for the 
modified equations that in particular satisfy the wave coordinate condition along 
the Cauchy hypcrsurface S; i.e., Q^|t=o = 0. As mentioned in the statement of 
Theorem 15. II these data launch a solution of both the modified equations and the 
Einstein equations. As we have discussed previously, this fact follows from the 
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fact that the condition = holds in 2?(S). We now formulate this result as a 
proposition. 



Proposition 5.6.1 (Preservation of the Wave Coordinate Condition). Let 

{Yi — {x £ M. \ t = 0}, g^i^, K_jf^,p,w'), {j,k = 1,2,3), be initial data for the Euler- 
Einstein system p. 2. lap - p.2.1cp that satisfy the constraints p.2.8ap - p.2.8bp . 
Let (g^,.|t=o,9tg^z.|t=o,p|t=o,'«''k=o), = 0, 1,2,3), he initial data for the modi- 
fied equations (|5.1.5a[) - (|5.1.5c[) that are constructed from the data for the Einstein 
equations as described in Section \5.3\ In particular, we recall that the construction 
of Section \5.3\ leads the fact that Q^\t=Q = where is defined in (|5.1.2b[) . Let 
{M, gfii,,p,uf^) be the maximal globally hyperbolic development of the data for the 
modified equations. Then = in 2?(S) = A4. 



Remark 5.6.1. To prove Proposition 15.6.11 we do not need the full assumption 
that the equation of state is of the form p = c^p; we only need the assumption 
p = p{p) and the assumptions described in Section [3] 



Proof. This is a rather standard result whose main ideas can be traced back to 
[CB52| . The sketch of a proof given in |RS09[ Proposition 5.6.1] can be easily 
adapted to apply to the system we are studying here. The basic idea is to show 
that for any solution to the modified equations, the quantities = — are 
solutions to a homogeneous system of wave equations with trivial initial data, i.e., 
that Qfi\t=o = dtQfj.\t=Q = 0. The fact that = tlren follows from the standard 
uniqueness theorem for such systems. □ 



5.7. Matrix-vector notation for the fluid variables. To streamline our no- 
tation, it will often be useful to write the fluid equations (j5.4.3d[) - (j5.4.3ep using 
PDE matrix-vector notation. To this end, we introduce the array of fluid variables 
W, the array of decay-inducing inhomogeneous terms b, and the array of error 
inhomogeneous terms b^, which arc defined by 





( P-P\ 




/ 


- 2)oju^ 

- 2)^1*2 




( ^\ 

A2 


(5.7.1) W = 


u 


, dcf 




, dcf 
DA = 




\ ) 






-2)wm3 j 







In the above formulas, x ~ 3c^, and A, A^ are defined in (j5.4.7dp - (|5.4.7ep . 
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We also introduce the 4x4 matrices A^, which are defined by 
(5.7.2a) 



(5.7.2b) 



A' 





nio 




p 1 


/ cl 




n2o 




p 1 


I 




n3o 




p J 



-(l + c2)PHl -(1 











' Mo V 



(1 



(l + c2)P 

{l+ci)P 



Vi 



l + c2)P 



n2i 
n3i 



c2)P 
1 












-c2)PH3\ 






and analogously for for A^, A? . For future use, we also calculate that the matrix 



A° satisfies det(A°) (u°)^| (m°)^(1 - c^) - c^5°° |, and its inverse is given by 
(5.7.3) 







p 








p 








p 



j20 
T30 



n 



Using this notation, the fluid equations (j5.4.3dp - (|5.4.3eP can be written in the 
following compact form: 



(5.7.4) Af^dfiW = \i + hc^. 

We remark that wc have split the inhomogcncous term into two pieces to facilitate 
our analysis in later sections. 

In Section 19.21 wc will provide estimates for the time derivatives of the fluid 
quantities. Therefore, as a preliminary step, we isolate them in the next corollary. 

Corollary 5.7.1. Let (P, u^, m^, w'^) he a solution to the modified relativistic En- 
ter equations (j5.4.3d[) - (j5.4.3cp . where is defined by (j5.4.4p . Then the time 
derivatives of the fluid quantities can be expressed as follows: 



(5.7.5a) 
(5.7.5b) 
(5.7.5c) 



dtu" = A'", 
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where 



(5.7.6a) 



X { - U^daP - (1 + C^)Paa^/'^ - ^i^J-^U^ffb^M^aau" 

(5.7.6b) 

A'" = i-(AO - _^i_n""A' - _^L_n"''a„p|, 

yOl (l + c2)F (l + c2)P i' 

(5.7.6c) 

A'^- = -i{ - A„^ + ....5°^ - ^_£L_n"^ A' - ^^^Ln'^.a^P - u^d.u^}. 

Proof. To obtain (|5.7.5ap . simply multiply both sides of (|5.7.4p by the first row 
of (AO)-i, use ((5Jl|) . (|5.7.2bp (plus the analo gous, but unwritten, formulas for 
j4^,j4'^) and (|5.7.3p . and consider the first component of the resulting expression. 
Equation (|5.7.5cp then follows from isolating dtu^ in equation (|5.4.3ep and using 
()5.7.5ap to substitute A' for dt{P — p)- Equation (|5.7.5bp similarly follows from 
([CT^ and (|5.7.5ap . 

□ 



6. Norms and Energies 

In this section, we define the Sobolev normf0 and energies that will play a central 
role in our global existence theorem of Section [TT] They are designed with equations 
(|5.4.3ap - (|5.4.3ep and Theorem [521 in mind. Let us make a few comments on these 
quantities. First, we remark that in Section [TUl we will show that if the norms are 
sufficiently small, then they are equivalent to the energies; i.e., the energies can 
be used to control Sobolev norms of solutions. The reason that we introduce the 
energies is that their time derivatives can be estimated via integration by parts/ 
the divergence theorem. Next, we recall that the background solution variables 
(P, u^, M^, u^) satisfy (P, m^, u^, u'^) = (p, 0, 0, 0) wherep > is the initial (constant) 
pressure, and the rescaled pressure P is defined in (|5.4.2p . The quantity Sp-p^u*;N 
introduced below in ()6.1.2b|) measures the difference of the perturbed fluid variables 
(P, M^, M^, u^) from the background {p, 0, 0, 0). We also follow Ringstrom |Rin08j by 
introducing scalings by e"^, where oc is a number, in the definitions of the norms 
and energies. The effect of these scalings is that in our proof of global existence, 
a convenient and viable bootstrap assumption to make for these quantities is that 
they are of size e, where e is sufficiently small. Finally, we remark that the small 
positive number q that appears in this section and throughout this article is defined 
in (|8.1.2p below, and we remind the reader that hjk == e^^^gjk, P == e^^^+^^^^p. 



Technically, Sp— p,u*;]v is a norm of the difference between the perturbed fluid variables and 
the background FLRW fluid solution (P,u^ ,u'^ ,u'^) = (p, 0,0,0). 
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6.1. Norms for g^^^, P — p, . 

Definition 6.1.1. Let be a positive integer. We define the norms Sg„„+i;M{t), 
Sgo.;Nit), Sh,,;Nit), Sg-Nit), UN-lit), Sp-p^u'-Nit), and S^it) as follows: 

(6.1.1a) 

Sg„„+uN =e^''\\dtgoo\\H-+e'^''\\goo + l\\H-+e^'-'^''\\d3oo\\H-, 
(6.1.1b) 

(6.1.1c) 

3 

Sh.,;N = {^'^WdMlH- + \\dh,k\\H--^+eS''-'^''\\dh,k\\H^ 



(6.1. Id) Sg-N — ' Sggg-^-l-N + Sgg^-N + Sh,,; 



(6.1.2a) C/^_/='e(l+')"(^||^.^■f^„-.)'^^ 
(6.1.2b) Sp.p.u-.-N e^^Wu^HN + \\P - pWh- , 

(6.1.3) Sn = Sg-N + Sp-p^u'iN + Un-1- 

Note that as discussed in Section 11.21 the norm Un-i in (|6.1.3p controls the 
lower-order derivatives of the with larger weights than the norm Sp^p^u* -,n ■ 

6.2. Energies for the metric g. 

6.2.1. The building block energy for g^,^. The energies for the metric components 
will be built from the quantities defined in the following lemma. They are designed 
with equations (|5.4.3ap - (|5.4.3cp in mind. 

Lemma 6.2.1. [Rin08| Lemma 15] Let v be a solution to the scalar equation 

(6.2.1) Dgv^ ocHdtv + ^H^v + F, 

where Og = g^'^Oxd^, a > and |3 > 0, and define 5(y dv] > by 

(6.2.2) 

,.[v, dv] i / { - g'"{dtv)^ + g''\dav){dbv) - 2yHg''%dtv + bH\^] d'x. 

Then there are constants!] > 0, C > 0, 5 > 0, andy > 0, withx] andC depending 
on (X, |3, Y and 6, such that 



(6.2.3) 

implies that 
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(6.2.4) £f ,.[v,dv]>C {dtvf+g''''{dav){dbv) + C^y^v^d^x, 

where Cf^y) = i/ y = and C(y) = 1 i/ y > 0. Furthermore, i/ |3 = 0, then 
J = b — 0. Finally, we have that 

(6.2.5) 

j^{Ely,,)[v,dv]) < -r]H£f^^,^[v,dv]+ J^^ { - {dtv + yHv)F + A£,^y^,)[v, dv]} d^'x, 
where 

(6.2.6) 

A£;(Y,5)b,at.] = -yHidag^'jvdbV - 2y H{dag°'')vdtv - 2yHg'"^{dav){dtv) 
- idag'"'){dtv)^ " {dag'"')idbv)idtv) - \idtg'°)idtv)^ 

ab I cr„ab Wq „,\/Q „ \ ^, tt t Pi „00\„. a „, ^,Tjf„00 



+ (^2^*5 + ^5 j idav){dbv) - yHidtgnvdtv - yH{g"^ + l){dtvY 

Proof. A proof based on a standard integration by parts argument (multiply both 
sides of equation (|6.2.ip by —{dtv +yHv) before integrating by parts over T^) is 
given in Lemma 15 of jRin08| . In particular, we quote the following identity: 

(6.2.7) 

^^{£f^^^^[v,dv]) = J^^{- {a-y)H{dtv)^ + (5 - ^-y<x)H\dtv - ^yH\^ 

- (1 + y)Hg''\dav) (dbv) - [dtV + yHv)F + A£.^^y^^)[v, dv] } d^'x. 

□ 

6.2.2. Energies for the components of g. 

Definition 6.2.1. Wc define the non-negative energies Eg^^+i-^M{t), Egg^.N{t), 
Eh„;N{t), and Eg.]y{t) as follows: 

(6.2.8a) 

Ko+i;N= E e2«^f(2.,„„,,„„)[as(.goo + l), 9(95300)], 

\a\<N 

(6.2.8b) 

[S\<N J = l 

(6.2.8c) 

3 

EL;N= E { E e''''4.o)[0'5aA^)] + J / csH'{dsh,,yd^, 

N^l^ i\r ^ 1 1 JT^ 



t\<N j,k=l 

(6.2.8d) Eg-N =* -Egoo+i;Ar + Eg^^-N + Eh,,-N, 
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where 



(6.2.9a) 
(6.2.9b) 
(6.2.9c) 



7 dcf _20 

hjk = e gjk, 

dcf 
CS ^ 0, 



(i,fc = 1,2,3), 
if |a| = 0, 
if |a| > 0, 



dcf 



1, 



and (yoo, 5oo), (Yo*, §0*), and (y,,,6„) = (0,0) are the constants generated by 
applying Lemma [6. 2. II to equations (|5.4.3a|) - (|5.4.3c|) respectively. 

In the next lemma, we provide a preliminary estimate of the time derivative of 
these energies. 

Lemma 6.2.2. jRS091 Lemma 6.2.2] Assume that g^^, {fijiy — 0,1,2,3), is a so- 
lution to the modified equations ()5.4.3a|) - (|5.4.3cp . and let Eg^g-^-i-N, Eg^^-^N, and 
Eh,,-N be as in Definition \ 6.2.l\ Let [Hg^ds] denote the commutator of the op- 
erators Hg and ds- Then under the assumptions of Lemma \6.2.1[ the following 
differential inequalities are satisfied, where A^.j^.^ ,5-) [•, 9(-)] is defined in (j6.2.6p . and 
the constants (yoo, §00), (Yo*, §0*), '^'^d (y**,5**) = (0,0) are defined in Definition 



(6.2.10a) ^(El^^^,.^^) < {2q - ^^^)HEI^^^.^^ + 2q{u - H)EI 



'soo + l:Af 




X 



{95A00 + [Ug, ds\{gm + 1)} 'f' 



X 





X { - 2Hds{g''''Tajb) + 9s Ao, + [fig, ds]gaj} d^x 




3 



+ E E / e2(^-i)^^A£;(^„.,5o.)[555o„5(<95ffo,)]d'a;, 
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(6.2.10c) 



3 

- J2 e^'^'\dtdsh,k + y*. H^sh,k}{^s^,k + [^g,^s]h,k}d• 
\s\<Nj,k=l•^'^'' 

3 

/t3 

|5|<Arj,fc=i-^^ 

1<|5|<W^T3 

Proof. Lemma 16.2.21 follows easily from definitions ()6.2.8ap - ()6.2.8d[) , and from 
(I6X5|) . □ 

The following corollary follows easily from Lemma 16.2.21 definitions (|6.2.8al) - 
(j6.2.8cp . and the Cauchy-Schwarz inequality for integrals. 

Corollary 6.2.3. |RS09| Corollary 6.2.3] Under the assumptions of Lemma \6.2.2\ 

we have that 



(6.2.11a) 

+ ^soo+i;W Yl e'''ll[ng,55](.9oo + 1)11^2 

\a\<N 

+ E e^''^i'^£;(Yoo,5oo)[55(5oo + l),c'(a5ffoo)]||Li, 

\a\<N 



3.2.11b) 

J^{EI,,n) < [2(<Z-l)-Tlo*]i^i?L;A^ + 2(9-l)(^-^)^L;^ 

3 3 
3 

+ Syo,;N E E^''"''"IIP5'^'S]50.IU= 

\a\<N j=l 
3 

+ E E^''''"''''ll^£;(Yo..6o.)[5s5o„a(a5ffo,)]IUi, 
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3.2.11c) 

d ^ 



dt 



+ Sh,.-N ^ ^ \\[^g.ds]hjk\\L'' 

\a\<N i,k=i 

3 

\a\<N j,k=l 

where the norms S'goo+i;7Vj 'S'go,;jVi Sfi,,-N defined in Definition \6.1.1\ 

6.3. The equations of variation and energies for P — p, uK In this section, we 
define energies that are useful for studying the Euler equations ()5.4.3d| - ()5.4.3ep . 
Now in order to estimate the spatial derivatives of the fluid variables, we will of 
course have to differentiate the equations (j5.4.3dp - (|5.4.3ep . or equivalently, equa- 
tion (|5.7.4p . The derivatives are solutions to the linearization of (|5.4.3dp - ()5.4.3ep 
around the background variables (P, m^, it^, u'^). We refer to the linearized system 
as the equations of variation, while the unknowns (P, u^, u^, u"^) are the variations. 
More specifically, the equations of variation in the unknowns (P,v} ,v? ^ir") corre- 
sponding to the background (P, u^, w^, u"^) are defined by 

(6.3.1a) u"5„P +{l + cl) (—) Puadtu" + (1 + cDPOaii" = ^, 

V Uq / 

(6.3.1b) + ff"a„P ={>(- 2)uju^ + 0\ 

(1 + c|)P 

where > is such that gapu^u^ — —1. In our applications below (P, ii^ , ii'^ , ii^) 
will be equal to (^ds{P ~p), dsu^ , dsu^, d^u^^ , where cS is a spatial derivative multi- 
index. The terms 3^, (>^— 2)aj'u^, & denote the inhomogeneous terms that arise from 
differentiating (j5.4.3dp - (j5.4.3ep . Note that we have split the inhomogeneous term 
in ()6.3.1b| into two pieces; the (>? — 2)uju^ piece is responsible for creating decay in 
the ii^ variable, while in our applications, & will be an error term. 

Using matrix-vector notation, we can abbreviate the equations ()6.3.1ap - (j6.3.1bp 

as 

(6.3.2) Af^dpW = I, 
where the A^" are defined in (|5.7.2ap - ()5.7.2bp . and 

(6.3.3) W"^' {P,u\u^,uY, 

(6.3.4) I = (i?,©i,0^©3)T^ 

To each variation (u^, u^, it'^), we associate a quantity -uP defined by 

(6.3.5) uO''^' ^Luau". 

uo 
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This quantity appears below in the expression (|6.3.13p . which defines our fluid 
energy current. In our analysis, we will need the following lemma, which essentially 
states that vP is a solution to the linearization of (|5.4.8p around the background 
(P,u1,m2 

Lemma 6.3.1. Assume that {P, ii^ , , ii^) is a solution to the equations of vari- 
ation (|6.3.1a|) - (|6.3.1bp corresponding to the background {P, , u"^ , u^) , and let 

vP —-^Uaii"' be as defined in (j6.3.5p . Then is a solution to the following 
equation: 



(6.3.6) 

where 



(l + c2)P 



(6.3.7a) 
(6.3.7b) 



(/i = 0,l,2,3). 



0U ^ _ 



Proof. By the definition of vP , the left-hand side of (|6.3.6|) is equal to 



(6.3.8) 



tio (1 + ci)P 



— 



On the other hand, contracting equation (j6.3.1b[) against Uj and using the iden- 
tity Mali"" = — uqII^", we conclude that 



(6.3.9) 



UaU°'da'U°' — Uq- 



-n°"daP = (>? - 2)ujUaU'' + Ua^". 



'{l + cf)P' 

Multiplying ()6.3.9p by — - and using (j6.3.8p . we arrive at (|6.3.6p . 



□ 



In the next lemma, we provide detailed information about the structure of the 
inhomogeneous terms appearing in the equations of variation and in Lemma l6.3.1l 
We again split the terms into two pieces to facilitate our analysis in later sections. 



Lemma 6.3.2. Let W = (P ^ p,u^ ,u'^ ,u^)"^ be a solution to the relativistic Euler 

equations (j5.7.4p . Then {P, , vP , ii'^)"^ {ds{P ~ p) , dsu^ , d^u^ , dsu^)^ is a so- 
lution to the equations of variation (j6.3.2p with inhomogeneous term I that can be 
expressed as follows: 



(6.3.10) 



I = dsh + bA5, 
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whe 



(6.3.11a) 

b =^ (0, - 2)u\ {k ~ 2)u^, {k - 2)u'^f, 
(6.3.11b) 

hAa = { A^^s [( A° ) " ^ b] - d^h} + dsliA^^^hA] 

+ A^^^iA^y^A^dadsW - ds[{A")-^A''daW]y 

Furthermore, iiP = —:^UaU°' is a solution to equation (|6.3.6p with inhomoge- 
neous term defined by 



(6.3.12) 25^ 



Proof. Equations (|6.3.10p - (|6.3.11bp are a straightforward decomposition of the 
inhomogeneous term I = A°ds{iA°)-^ A'^df.'W} + A°[{A°)-^At'd^,d3]W 
= A°ds{{A°)-^{h + h A)} + A°[(A°)-M°9a,9a]W where [•,•] denotes the com- 
mutator. The relation (j6.3.12p follows directly from (j6.3.7b|) . 

□ 



6.3.1. The fluid energy currents. To each variation W = {P,u ,u ,ir) , we asso- 
ciate the following energy current, where vP —-i^UaU'^ '■ 



(6.3.13) > (TT^'"' + + ^^if^5./.u"u^. 

Currents of the form (|6.3.13p are the building blocks for some of our fluid energies 
(see (|6.3.15p ). We remark that similar currents were used in |Chr07bj . |Spe09a| , 
|Spe09b| . 



Remark 6.3.1. We sometimes write J^[W,W] to emphasize that depends 
quadratically on the variations. 

6.3.2. The divergence of the fluid energy current. Let W ~ {P, ii^ , ii^ , ii^)^ be a 
solution to the equations of variation (|6.3.1ap - (j6.3.1bp . Then using the equations 
(|6.3.1ap - (j6.3.1b[) and (|6.3.6p to replace derivatives of W with the inhomogeneous 
terms, an omitted, tedious computation gives that the divergence of J'^ can be 
expressed as follows: 
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(6.3.14) 

9^(>[W,W]) 



dt 



(1+C2)P 

(1 + c2)FwO 



da 



p2 



(1+C2)P. 



da 



(1 + c^)Pm° 
(1 + c?)Pw" 



{dagmW) 



.goo(w°)2 + 2.goaM°u" 



{dagira)u U 



I •, rn 



{l + cl)Pu^ 



9abU U 

-{dagob)u°u'' 



2(l + c2)PwO 



+ 



2(U 



'Ua 
.Uq 



idtgoa)u°u'' -r - 
{l + cl)P 



(l+c2)P(u"-l) 



idtgabKu' 



(X - l)LOgabU°'u' + 



{dtgab - 2ujgab)u''u^ 
2^ 



[l + cl)P 



p 



<0 



H 5 gm'o u 



A{\ + cl)P 



A'^o<e t/iai </ie right-hand side of (j6.3.14p lioes not depend on the derivatives of W; 
this property is essential for closing our energy estimates for the top derivatives of 
the fluid variables. We remark that we have organized the terms on the terms on 
the right-hand side of (|6.3.14p in a way that will be useful for proving inequality 
(111.1.311 . 

6.3.3. The definition of Ep-p^u*;N ■ In our analysis, we will make use of several 
different norms and energies for the fluid. The energy Ep^p^u*;N will control all 
spatial derivatives of all of W (P—p, v} ^v?, u^)^, while the norm f/jv-i defined in 
(|6.1.2ap will control the lower derivatives of with larger weights. These 

larger weights lead to better L°° decay for the lower-order derivatives of {u^ , u^) 
than for the top order derivatives; this improved decay plays an essential role in 
our analysis. We now proceed to the definition of the non-negative (see inequality 
(|10.0.46f|) ) quantity Ep-p,u-:N- 

Definition 6.3.1. Let be a positive integer, and let W =* [P — p,u^ ,v? ,vr')'^ 
be the array of fluid variables. We define the fluid energy Ep^p,u*:N{t) > by 



(6.3.15) 



^P-p,u':N 



V / j''[d3W,dsW]d''x, 



|q|<JV 



where j°[ds'W,dsW] is defined in (l6.3.13|) . 

In the next corollary, we provide a preliminary estimate for the time derivatives 
of Un-1 and Ep^p^^'-.N- 
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Corollary 6.3.3. Let W == (P — p,u^ ,u^)"^ be a solution to the relativistic 
Euler equations ()5.7.4p . Let U^-iit) and Ep^p^u*:N{i) be the fluid norm and en- 
ergy defined in (|6.1.2aP and (|6.3.15p respectively. Then the following differential 
inequalities are satisfied: 

(6.3.16) 

< 2 (>. - 1 + g) ^e(i+')^^f/^_i + 2e(i+«)^^C/^_i ^ l|A-||««-., 

(6.3.17) < ^ ||a^(>[9„-w,a5W])||^,. 

|q|<W 

Proof. To prove (|6.3.16p . we use the definition (|6.1.2ap of Um-i and equation 
(|5.7.5cp (differentiated with Sg) to conclude that 

(6.3.18) 

^(C/2,_,) =2(l + g)^[/2,_,+2e2(i+9)f^ E E / {dsn-){{K - 2)u:dr.u- + dr.^'= 

Inequahty (|6.3.16p now fohows from (|6.3.18p . the definition of Un-i, and the 
Cauchy-Schwarz inequahty for integrals. 

To prove (|6.3.17p . we use the definition (j6.3.15p of Ep ^ and the divergence 
theorem to conclude that 

(6.3.19) 1(£;2__„.^^) ^ ^ f dt{j''[dsW,dsW])d^x 

= E / d^{j''[dsW,dsW])d^x, 

from which (I6.3.17P easily follows. 

□ 

6.4. The total energy Eat. 

Definition 6.4.1. Using definitions (|6.1.2ap . (|6.2.8dp . and (|6.3.15p . we define Eat, 
the total energy associated to g, P, u, as follows: 

(6.4.1) Eat = Eg-N + Ep^p^u';N + Un-1- 

7. Linear- Algebraic Estimates of and g^'', (/i, = 0, 1, 2, 3) 

In this section, we recall some linear-algebraic estimates of g^i, and that were 
proved by Ringstrom. In addition to providing some rough L°° estimates that we 
will use in Sections [S] and 1101 the lemmas will guarantee that g^^, is a Lorentzian 
metric. This latter fact has already been used in our statement of the conclusions 
of Theorem 15.11 
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Lemma 7.0.1. |Rin08[ Lemmas 1 and 2] Let gf^^, be a symmetric 4x4 matrix of real 
numbers. Let {g\,)jk be the 3x3 matrix defined by {gi,)jk ~ gjk, and let {g^^Y'' be 
the 3x3 inverse of {g\,)jk- Assume that goo < and that {g\,)jk is positive definite. 
Then g^^ is a Lorentzian metric with inverse g^'^ , g"'' < 0, and the 3x3 matrix 

{g'^y^ defined by {g'^y^ = g-'^ is positive definite. Furthermore, the following 
relations hold: 



.goo - d"^ ' 



(7.0.1a) 3™ = 

(7.0.1b) 

.900 - " 

(7.0.1c) = —^{g-^y^g,,, (j = 1, 2, 3), 



(P - goo 



where 



(7.0.2) 



ig\, ^T'^goagob- 



□ 

The estimates in the next lemma are based on the following rough assumptions, 
which we will upgrade during our global existence argument. 

Rough Bootstrap Assumptions for g^y : 

We assume that there are constants t] > and ci > 1 such that 

(7.0.3a) |ffoo + l|<Tl, 

(7.0.3b) c^^SabX^X^ < er'^^^gabX^'X'' < ci6abX''X'', y{X\X^, X^) G M^ 



(7.0.3c) 



For our global existence argument, we will assume that r| = r\mim where rimin is 
defined in Section [01 



Lemma 7.0.2. [RinOSj Lemma 7] Let g^^ be a symmetric 4x4 matrix of real 
numbers satisfying (|7.0.3ap - (|7.0.3cp . where > and < g < 1. Then g is a 
Lorentzian metric, and there exists a constant r\o > such that < Tl t]o implies 
that the following estimates hold for the its inverse g^'^ : 



(7.0.4a) 
(7.0.4b) 

(7.0.4c) 



I<?™ + 1| <4ri, 



-2n 



\ a=l 



\ a=l 



\g"'goa\<2c^e-'''Y.\9oa\^ 



(7.0.4d) -^S''''XaXi, < e'^g'^'^XaXk < ^S'^'XaXt, V(Xi, X2, X3) e M.\ 



2 

3ci 
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□ 

8. The Bootstrap Assumption for Sn and the Definition of Af,rimm, 

AND q 

For the remainder of the article, N denotes a fixed integer subject to the require- 
ment 

(8.0.5) 

N > 3 (This is large enough for all of our results except some of the conclusions of Theorem 112. ip , 
(8.0.6) 

N > 5 (This is large enough for the full results of Theorem 112.11 to be valid) . 

We require TV to be of this size to ensure that various Sobolev embedding results 
and the conclusions of the propositions in Appendix [A] are valid. 

In our global existence argument, we will make the following bootstrap assump- 
tion: 

(8.0.7) Sn < e, 

where S^v is defined in (|6.1.3p . and e is a sufficiently small positive number. Observe 
that Sat measures how much {g,p,u) differs from the FLRW background solution 
(]j,p,u) derived in Section |4l In particular, S^r = for the background solution. 

8.1. The definitions of r|,„i„ and q. 

Definition 8.1.1. Let us apply Lemma [6.2.11 to each of the equations (|5.4.3ap - 
(|5.4.3cp . denoting the constant r| produced by the lemma in each case by riocTio*, 
and T)** respectively. Furthermore, let r)o be the constant from Lemma [7.0.21 We 
now define the positive quantities (recalling that < x < 1 when < Cg < -^/l/B) 
r[min and q by 

(8.1.1) r],nin imin{l,Tio,Tioo,'no*,'n*,}, 

(8.1.2) 9 =' ^min{r|„ij„, 1 - 

The constants r\min and q have been chosen to be small enough to close the boot- 
strap argument for global existence given in Section lll.21 In particular, inequality 
(|7.0.4a|) . with ri < rimi„, guarantees that the energies £(^^5) [■, 9(-)] for solutions to 
(|5.4.3a|) - (|5.4.3c|) have the coercive property (|6.2.4p . 

Remark 8.1.1. By Sobolev embedding and Lemma 17.0.21 if e is small enough, 
then the assumption S^r < e implies that there exists a constant C > such that 

(8.1.3a) |5oo + l|<Ce, 

(8.1.3b) j2\9oj\<CeeM-''^''. 

i=i 

Therefore, if e is sufficiently small, the inequalities (|7.0.3a|) and ()7.0.3c|) are an 
automatic consequence of (|7.0.3b|) and the assumption Sn < £■ 
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9. SoBOLEv Estimates 

In this section, we use the bootstrap assumptions of Sections [7] and [8] to deduce 
estimates of g^^, g'"^, P, li^, and the nonhnear terms in the modified equations 
(|5.4.3a|) - ()5.4.3e|) , (|5.4.8p . The main goal is to show that the error terms are smaU 
compared to the principal terms, which is the main step in closing the bootstrap 
argument in the proof of Theorem 111.11 The primary tools for estimating these 
terms are standard Sobolev-Moser estimates, which we have collected together in 
the Appendix for convenience. 

9.1. Estimates of g^^, g'^'^. In this section, we state the first proposition that 
win be used to deduce the energy inequalities of Section fIl.il The proposition was 
essentially proved in |RS09[ Proposition 9.1.1], using the ideas of |Rin08| Lemmas 
9,11,18,20]. We don't bother to repeat the proof here, since similar arguments are 
used below in the proof of Proposition 19.2.11 We remark that the proof makes use 
of the lemmas stated in Section [71 

Proposition 9.1.1. |RS09i Proposition 9.1.1] Let N > 3 be an integer, and assume 
that the bootstrap assumptions (|7.0.3ap - (|7.0.3cP hold on the spacetime slab [0, T) x 
T'^ for some constant Ci > 1 and for r\ = T)mm- Then there exists a constant e' > 
and a constant C > 0, where C depends on N, ci, andr\mim such that if 8^(1:) < e' 
on [0,T), then the following estimates also hold on [0,r), where hjk = e^'^^gjk ■ 



(9.1.1a) ILgoolU- < 2, 

(9.1.1b) llffo.lk- < Ce(i-'^)"5g;^, 

(9.1.1c) Wgjkh^ <Ce^^, 

(9.1.2a) < Ce-?"5<,;Ar, 

(9.1.2b) ||3°1l~<5, 

(9.1.2c) ||/° + l||^« <Ce-?"S<,.^, 

(9.1.2d) Wg"" + <Ce-''''Sg,N, 

(9.1.2e) <Ce-2", 

(9.1.2f) <Ce-^"Sg.N, 

(9.1.2g) \\d3'''\\L^<Ce-^"Sg,N, 

(9.1.2h) < Ce-(i+«)^^5<,;^, 

(9.1.2i) < Ce-(i+')^^5,;Ar, 

(9.1.3a) \\dtg,k - ^ug^kWH- < Ce^^-''^''Sh,.;N, 

(9.1.3b) Wdtgjk - 'i^gjkWcl ^ Ce^^~''''"Sh„;N, 

(9.1.3c) \\dtgjk\\cl<Ce^", 
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(9.1.4a) Wg^'dtgak - 2w5i||ff« < Cer'^^'Sg-.N, 

(9.1.4b) Wdtgak - 2uj5i\\L^ < Ce~'^''Sg,N, 

(9.1.5a) + 2u;g^'=||^« < Ce-^^+'^^''Sg,N, 

(9.1.5b) 119(5^"'= + 2ojg^''\\L^ < Ce-^^+'i^'' Sg,N , 

(9.1.5c) ||9t5™IU- <Ce-«"5<,;Ar, 

(9.1.5d) <Ce-(i+')"5,;w, 

(9.1.5e) \\dtg'''\\L^<Ce-^'\ 

(9.1.6a) ||5'"r,,fc||H« < Ce(i-«)^^^^„;jv, 

(9.1.6b) h'^'TaAn--^ <CSh..-N. 



In the above estimates, the norms Sh -n o,nd Sg-M are defined in Definition 

[223 

□ 



9.2. Estimates of P, u^, and the error terms. In this section, wc state and 
prove the second proposition that wiU be used to deduce the energy incquahties of 
Section nm 

Proposition 9.2.1. Let N > 3 be an integer, and let [gf^^, P, u^), (/i, = 0, 1, 2, 3), 
be a solution to the reduced equations (|5.4.3aP - (|5.4.3eP on the spacetime slab 
[0,T) X T'^. Assume that the bootstrap assumptions (|7.0.3ap - (j7.0.3c|) hold on the 
same slab for some constant Ci > 1 and for r\ = r\min- Then there exists a constant 
e" > and a constant C > 0, where C depends on N,ci, and r\mini such that 
if SN{t) < e" on [0,T), then the following estimates also hold on [0,T) for the 
quantities A^.^i,, AcoOi 'I'^d Ac,oj, defined in (|5.4.17aP - (|5.4.17cP and (|5.4.19aP 
- (|5.4.19bl) .- 



(9.2.1a) 


||Aa,oo 


HN 


< Ce-^^^Sl^, 


(9.2.1b) 


I|Aa,oj 


HN 




(9.2.1c) 


l|AA,,fc 


HN 




(9.2. Id) 


l|Ac,oo 


HN 




(9.2. le) 


l|Ac,oj 


HN 
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For the fluid quantities, we have the following estimates on [0,r) : 

(9.2.2a) llPlUoo < C, 

(9.2.2b) - lll^iv < Cc-'^'^Sn, 

(9.2.2c) ||u"-1||loo < Ce-«"Sw, 

(9.2.2d) ||uo + 1||h" < Ce-«^Sw, 

(9.2.2e) ||«o + 1||l- < Ce-'^^'SAT, 

(9.2.2f) ||uj||^iv-i < Ce(i-9)"SA,, 

(9.2.2g) \\u,\\l^ < Ce(i-«)"SA., 

(9.2.2h) < Ce^Sjv. 

For the time derivatives of the fluid quantities, we have the following estimates 
on [0,T) : 

(9.2.3a) ||atP||H"-i < Ce-'i^SN, 

(9.2.3b) WdtPh^ < Ce-«"SAr, 

(9.2.3c) ||9tu°||^iv-i < Ce-^^Sw, 

(9.2.3d) ||9tw°||L- <Ce-«"SAr, 



(9.2.3e) \\dtu^\HN-i <Ce-(i+«)"S; 

(9.2.3f) llStu^'IU- <Ce-(i+«'"S^ 



(9.2.3g) \\dtUo\\H«-^ < Ce-^^SAT, 

(9.2.3h) Wdtuoh^ <Ce-'^''SN, 

(9.2.3i) Ptujllffiv-i < Ce(i-«)"S7v, 

(9.2.3j) llStUjIlioo <Ce(i-9)"SAr. 

For t/ie quantities A^^, defined in (j5.4.7ap - (j5.4.7cp . we have the following esti- 
mates on [0, T) : 

(9.2.4a) ||Aoo||h« < Ce-^^^^Sjv, 

(9.2.4b) I|Aoj||h« < Ce(i-29)^^Sw, 

(9.2.4c) \\A,k\H- < Ce-^'''SN. 

For the commutator terms from Corollarv \6.2.3\. we have the following estimates 
on [0,T) : 

(9.2.5a) \\[ng,ds]{goo + 1)||l^ < Ce-^^^^Sjv, 

(9.2.5b) \\[Dg,ds]goj\\L- < Ceti-^^jOg^^ 

(9.2.5c) \\[Dg,ds]hjk\\L- < Ce-29^^Sjv. 
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For the terms from Corollary \ 6.2.3\ where A^.^^ ,5)[t;, (9u] is defined in (|6.2.6p . 
we have the following estimates on [0, T) : 



e2?"||A^,(^„„,5oo)[9a(.9oo + 1), 9(35.900)] IIli < Ce-'"5<,„„+i;jvSjv, 



(9.2.6a) 
(9.2.6b) 
(9.2.6c) 

For the Christoffel symbol error terms defined in (|5.4.21ap - (j5.4.21f|) . we have 
the following estimates on [0, T) : 



e2(^-i)^||A£^(^„.,5o.)[55(<?o,),5(c^55o,)]||Li < Ce-'"^,„.;^S^, 
e"^'lA£.(o,o)[0,a(55/i,fe)]||Li < Ce-'>''Sh,.;NSN- 



(9.2.7a) 
(9.2.7b) 
(9.2.7c) 
(9.2. 7d) 
(9.2.7e) 
(9.2.7f) 



HA " I 
1^0 ol 

IIA " , 
W^j ol 

IIAo^oi 

l|Ao\., 

IIAA, 



,\\hn <Ce-^'+''^''Sg,N, 
J|h« <Ce-'^"Sg,N, 

For the error terms A, A^, A', and A'^' defined in (|5.4.7dp . (|5.4.7ep . (|5.4.9p . 
(|5.7.6ap . (j5.7.6bp . (|5.7.6cP we have the following estimates on [0,T) : 



\\A\\hn <Ce-^''SN, 



(9.2.8a) 
(9.2.8b) 
(9.2.8c) 
(9.2.8d) 
(9.2.8e) 

(9.2.8f) \\A'^\hn-i < Ce-(i+«)f^Siv. 

For the norm of the variation vP defined in ()6.3.5p . we have the following 
estimate on [0, T) : 



^7 



(9.2.9) 



For t/ie norms of the inhomogeneous terms ^s, &'t defined in (j6.3.11bp . 
()6.3.12p . we have the following estimates on [0,T) : 



(9.2.10a) 



(9.2.10b) 



el 

a 

el 

\4j 



( 



< CS 



N 



L2 



g-(l+g)n 

g-(i+9)n J 



(0 < |a| < N), 



(0 < |a| < N). 



In the above estimates, the 
defined in Definition \6.1.1\ 
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Proof. 

Remark 9.2.1. Throughout the remaining proofs in this article, we wiU make 
use of the resuhs of Lemma 14.1.11 the definitions of the norms from Section |6l 
the definitions (|8.1.ip . (|8.1.2p of r\,nin and q, and the Sobolev embedding result 
jjM+2(j3^ ^ C^'{T^)^ {M > 0). We will also use the assumption that 
Sn, which is defined in (|6.1.3p . is sufficiently small without explicitly 
mentioning it every time. Furthermore, the smallness is adjusted as 
necessary at each step in the proof. For brevity, we don't explicitly estimate 
how small S^r must be. We also remark that as discussed in Section [2.51 the 
constants c, C, C* that appear throughout the article can be chosen uniformly 
(however, they may depend on N) as long as Sat is sufficiently small. Finally, we 
prove statements in logical order, rather than the order in which they are stated in 
the proposition. 

Before beginning the proof, we observe the following updated version of the 
Counting Principle from (RS09| , which provides useful heuristic guidelines for many 
of the estimates. This tool is only intended to help guide the reader through the 
estimates; we provide complete proofs of many of the estimates. 



Counting Principle 



Consider a product which contains as factors metric components g^^^^ inverse 
metric components g^'^ ^ the first derivatives of these quantities, and the four-velocity 
it^. IfU denotes the total number of upstairs spatial indices among these factors, 
and D denotes the total number of downstairs spatial metric indices, then the 
expected contribution to the rate of growth/decay of the norm of the product 
coming from these terms is no larger than (i.e., we expect these terms to 

contribute at least this much decay). For purposes of counting, a spatial derivative 
dj of a metric or inverse metric component is considered to be a downstairs spa- 
tial index, while time derivatives of these quantities are neutral. We remark that 
by these criteria, hjk e^^^gjk, {j,k = 1,2,3), is an order 1 term that doesn't 
contribute to the decay rate. Furthermore, each factor in a product, excluding hjk 
but including dihjk, that is equal to one of quantities under an norm in defi- 
nitions (|6.1.1a|) - (|6.1.2ap contributes an additional decay factor of e~''^. Finally, 
the rescaled pressure P counts as an order 1 term. 

In this way, many of the estimates proved below in detail can be 
ascertained by counting spatial indices and making sure that at least 
one factor in a product contributes an additional decay factor of e 



Proofs of (|9.2.1ap - (|9.2.1eP : To prove (|9.2.1ap . we first recall equation (|5.4.17ap : 
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(9.2.11) 



A 



A, 00 



{g°y{{dt9oof - (Fooo)'} +5°°.9''°{2(ft3oo)(9tffoa + ^affoo) - 4rooorooa} 

+ 5™5'''{(5t5oa) (9*506) + (aa5oo)(a65oo) - 2rooaroob} 

+ g"''g"''[2{dtgoo){da9ob) + 2{dtgob){dagoo) ~ 2roooraob - 2roo6rooa} 

+ g'''>g"'[2{dtgoa){digob) + 2{dbgoo){dagoi) - 4rooar,o6} 

+ g''''g'"'{dagoi){dbgo.n) + ^5""( g'^'dtgai - 2c^Jf ){dbgom + d^gob) 

2 „ ' 



-g g {dagol + 9;.9oa)(9hgom + 9m5o6) 

1 



-( rt5a/-2c.<5f )( g'^dtgbrn~2u;5i, ). 

4 V ^ ' V ^ ' 

e^"g'^i'dth^i-2i^g"i'goi e'^'^ g'"- dtMrr^-2ugO' gob 

We now use Proposition IA-51 the definition (|6.1.3p of Sjv, Sobolev embedding, 
(jgX^ (19.1. 2dl) . mrM . (19X21 . (I9.1.2hl) . (19X21)) . (I9.1.3al) . (|gXBa . (l9.1.4aL 
(j9.1.4bp . and the relation F^q,!^ = \{d^j,gau + d^g^fi — dag^j,^) to conclude that 

(9.2.12) ||A^,oo||h" < Ce-2«^^52^. 

This proves (|9.2.1ap . The proofs of (j9.2.1bp - (|9.2.1eP are similar, and we omit the 
details. We also remark that slight variations of these inequalities were proved in 
|Rin08| Lemma 12]. 



Proofs of (|9.2.2ap - (j9.2.2hp : Inequality (|9.2.2ap follows from the definition ((03)) 

of Stv and the Sobolev embedding result \\P — p||l=° < C\\P — p||ijN < CSm- 
To prove (|9.2.2bp . we first recall equation ([5.4.4^ : 



(9.2.13) 



500 



'1 



goo 



as 2 



gabWu'' _ / goo + 1 
.900 ^ .900 



We now apply Corollary I A-41 Proposition I A-5[ and Sobolev embedding to conclude 
that 



(9.2.14) 

ll„,0 1 



\hn<C 
< C 



goaU" 



goo 



goaU" 
goo 



gabU"u'> goo + 1 



.900 



.900 



1 








1 


2 






+ 




+ 


+ 




' 1 


500 




^.900^ 




500 




^500^ 


HN-l J 



(9.2.15) 



+ (||5a&||L- + \\d3ab\\H'^-^)\K\\L'^\\du''\\HN-i + ||500 + 
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Using Corollary IA-31 (|9.1.1cp . the definition of Sat, Sobolev embedding, and in 
particular the estimate Hm^Hl"' < Ce^^^+'^-'^^Sjv, it follows that 

(9.2.16) \\u° - 1\\hn < CSjve-"". 

This proves (|9.2.2bp . Inequality (|9.2.2cp then follows from (j9.2.2b|) and Sobolev 
embedding. 

Inequalities (|9.2.2dp - (|9.2.2ep can be proved similarly using (j9.2.2bp - ()9.2.2c|) 
and the triangle inequality estimate 

(9.2.17) 

\\uo + = Wgoc^u" + < !l(.goo + 1)u"||h« + \\u° - + WgoaU^WnN. 

Inequalities (|9.2.2f|) - (|9.2.2hp can also be proved similarly using the relation Uj — 

gjaU°'. 

Proofs of (|9.2.4ap - (|9.2.4cp : To prove (|9.2.4ap . we first use equation (|5.4.7ap and 
Proposition I A-51 to arrive at the following estimate: 

(9.2.18) 

IIAoolk" < c{i|A^^oo!lH« + ||Ac,oo!Ih« + e-'^'+''''' 11.900 + 1\\h- 

+ ||P||l~||uo + 1||l-||wo - + ll^o + l|k-=||f'||L-||wo - l||L~||9wo||ff"-i 
+ ho + lliio - ||aP||H«-i + e-3(i+^')^^||P||L^ 11.900 + 

+ e-3(l+'=^)^^||goo + l||L=e|[5P||H«-i 

+ \uj- H\\\dtgoo\\H'^ + \uj~ H\\\goo + 

We now use (|9.1.2cp . (|9.1.2dp . (|9.2.1ap . (|9.2.1dp . (|9.2.2ap . (|9.2.2dp . (|9.2.2ep . the 
definition (|6.1.3p of Sjv, and Sobolev embedding to conclude that 

(9.2.19) ||Aoo||h« < Ce-29f^Sjv, 

which proves (|9.2.4ap . Inequalities (j9.2.4bp and ()9.2.4cp can be proved using similar 
reasoning; we omit the details. 

Proofs of (|9.2.6ap - (|9.2.6cp : To begin, we first recall equation (|6.2.6p : 
(9.2.20) 

As,^^^S)[v,dv] = -yH{dag'"')vdbV - 2y H{dag°'')vdtv - 2y Hg'"^{dav){dtv) 

~ {dag'^'^rnvf - {dag'''){d,v){dtv) - i(9,500)(9,^;)2 

+ [\dtg'''' + ^5"') (9a«)(a6z;) + (H - uj)g''\dav){dbv) 
-yH{dtg°'')vdtv - yi?(g°° + l){dtvf. 
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We now claim that the foUowing inequahty holds for any function v for which 
the right-hand side is finite: 

(9.2.21) 



where C(y) is defined in (|6.2.4p To obtain (|9.2.2ip . we use the Cauchy-Schwarz 
inequality for integrals, (|9.1.2dp . (|9.1.2e|) . p.l.2gp , (|9X2i)) . (|9.1.5bp . and (|9.1.5cp . 
Inequalities (|9.2.6ap - (|9.2.6cp now easily follow from definitions (|6.1.1ap - (|6.1.3p 
and (|9.2.2ip . 

Proofs of (|9.2.7ap - (|9.2.7f|) : To estimate Aq^o, we first recall equation (|5.4.21ap : 
(9.2.22) AgOo = ^/°9t5oo + 5°"5tgoa - ^ff^^^agoo. 



Using Proposition IA-5[ the definition (|6.1.3p of Sa?, Sobolev embedding, ()9.1.2bp . 
(|9.1.2cp . (j9.1.2hp . and (|9.1.2ip . it follows that 



(9.2.23) 



which proves (|9.2.7ap . 

The estimates (|9.2.7bp - (|9.2.7fp can be proved similarly using Propositions 19. 1 . T] 
and lA-51 



Proofs of (|9.2.8ap - (|9.2.8f!) : To prove (|9.2.8ap . we first apply Proposition IA-51 to 
equation (j5.4.7dp . concluding that: 



(9.2.24) 

\\h« 



A\\h« < c|||P|U..HA„"o|U^||5i.°||h"-i|| + \\P\\l^\\u^\\l^\\A^\\\h- 

A„"ol|L~|k°||L~||aP||^«-. + \\PU^\\A^\h^\K\\H« 
" ^L-\\A^\\\H- + \\A^\\\L-\K\\L-\\dP\\H«-^ 



\\P\\l 

1 
1 



J\uYL^\\dt9oo\\H- 



1 

Uo 







1 


+ 






HN-1 


Uo 



||9t<?0a|U~|k1|L^||a^.°|U« 
+ !|9t.90a||L-|k"||Lo=||u1|L=e a(- 



Wdtgooh^Wu'h^Wdu^WH'^^^ 

\\dt90a\\L^\\u°\\L^\\u^H- 

\\u"\\L'^\\u''\\L'^\\dtgoa\\H« 
Jdtgab\\L^\K\\L^\\u''\\H^ 



1 

Uo 
+ 



1 

Wo 



JK\\L^\\u'^\\L^\mdt9ab)\\H--^ 



-Wu^L^Wu'WL^Wdtgabh^ 
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Using Corollary IA-31 to estimate d(^':^ 
(|9.2.2ep . (|9.2.7ap - (|9.2.7f|) the definition 



(19. I.Hall . (19.1. acll . (I9.2.2bll - 

3]) of S^r, Sobolev embedding, and 



in particular the estimate < C||m°||^n-i < Ce-(i+«)"SAr, it follows that 



(9.2.25) 



\A\\hn < Ce-'^^S 



N- 



This proves (|9.2.8ap . The proofs of (|9.2.8bp - (|9.2.8f|) are similar, and we omit the 
details. 



Proofs of (|9.2.3al) - ( [9X311 ): 

Inequality (|9.2.3a[) . follows trivially from equation (j5.7.5ap and from (j9.2.8dp . 

(|9.2.3dP then follows from (|9.2.3ap and Sobolev embedding. The proofs of (|9.2.3cp 

- (j9.2.3f|) are similar. 

The estimates ( |9.2.3gP - ( |9.2.3j[ ) follow from the identity 



(9.2.26) 



(Ai = 0,l,2,3), 



Proposition IA-5[ the definition (|6.1.3|) of Sjv, Sobolev embedding, the estimates of 
Proposition and (|9.2.3cl) - (|9.2.3fl) . 

Proof of dMill): To prove (|9X9l) . we first recall equation (|9.2.27|) : 



(9.2.27) 



1 

Uq 



■ dcf 

U = UaU 



Using (19.2. 2e|) and ( |9.2.2gp , it follows that 



(9.2.28) |l?i°|U2 < 

which proves (|9.2.9p . 



Proo/s o/ (|9.2.10ap - (|9.2.10bp : To prove (|9.2.10ap . we first use equation (|6.3.11bp 
to deduce that 



(9.2.29) 



©I 



< \\A°d4{A'')-'h]-dsh\ 



|A"9s[(A°)-ib 



A 



A°{(^°)-iA°9„95W - ds[{A°)-'A''daW]} 



lL2 



L2 



See the remarks at the end of Section 12.41 concerning our use of notation for the 
norms of array- valued functions. 

We will estimate each of the three terms on the right-hand side of (|9.2.29p using 
the following estimates for W =^ (P — p, ,u^, u^)^ , 

b= (0,(>f-2)ui,(>f-2)u2,(><-2)w3)^,bA - (A, Ai, A2, A^)^, A^, and {A°)-^ : 
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(9.2.30) ||b||^« 



-1 < CS 



N 



I 

g-(l+g)n 
g-(l+g)n 
(l+g)n j 



V' 



( \ 



N 



(9.2.31) P°||l-<C 
(9.2.32) 



1 e 



(i-9)^Sa, e(i-'?)^^Sjv e(i-«)f^SA.\ 



/ 1 e(i-«)"Sjv 
e-(i+9)"Sw 1 



1 



(9.2.33) ||9(A' 



(9.2.34) 



0\-l I 



e ----ajv 



-1 < CS 



N 



e 



1 



-qU 



lib 



< CS 



N 



^-{i+q)n ' 

£-(1+9)^^ 

yg-(n-9)ny 



(9.2.35) 



\\W\\hn < CSn 
/e-(i+g)ns 



■>N 



1 



(9.2.36) IIAiIlioc, <C 



V 



-20 
-20 



-20 



e 







-(1+9)^8^ 




(9.2.37) 



l|5A'l 



if" 



-1 < CS 



TV 





-n 






e~ 






2n 








\ 







and analogously for A^,A^. All of the above estimates follow from repeated ap- 
plications of Corollary IA-3[ Proposition IA-51 the definition (|6.1.3p of Sjv, Sobolev 
embedding, the estimates of Proposition 19.1.11 (|9.2.2ap - ( |9.2.3jp , and (|9.2.8ap - 

To estimate the first term on the right-hand side of (j9.2.29p , we use Propositions 
IA-5l and lA-6[ together with the above estimates and Sobolev embedding, to conclude 
that for < lal < iV 
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(9.2.38) 



L2 



< cs 



e-(i+9)f^ ' 
e-(i+9)iJ 
yg-(i+9)ny 



where we write * to emphasize that we are performing matrix muhiphcation on the 
matrices of norms. 

We estimate the second and third terms on the right-hand side of ()9.2.29p using 
similar reasoning, which results in the following bounds: 



(9.2.39) 

||A°{(A0)-M'^aa9sW - 95[(A°)-iA"9,W] } 

< cpO|Uo. * {||(A°)-i|u=. * m^H--. + ||a(A")-i^«-. * UA'^Wh--^ 



< CS 



N 



-212 
-2f2 



(9.2.40) 



< ||A"||l,c * <; ||(A0)-1||lc. * ||bA||H« + ||9(A")-^||Hiv-i * llbAlU 



0\-l 



< CS 



N 



-{i+q)n 



Finally, adding (|9.2.38p - (|9.2.40p implies that 
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(9.2.41) 



©I 

a 
©2 

\4j 



< cs 



N 



L2 



\^g-(i+9)ny 



which proves (|9.2.10ap . 

To prove (|9.2.10b[) . we use equation (j6.3.12p . the definition (|6.1.3p of S^r, Sobolev 
embedding, and in particular the estimates ^ Ce~'-^^'^^^SN,\\uj\\Lc^ < 

Ce(i-«)^^Siv, (|9.2.3dl) . ( [9T3j1 ), and (19.2. lOal) to conclude that 



< 






1 








uo 



L2 



< llu^llL-lia.Mal 



II VUO^ 



+ 



1 

Wo 



\a\<N 



\U Whn 



\Ua\\L° 



L2 



|5|<Ar 



This completes the proof of (j9.2.10bp . 

Proofs of (|9.2.5ap - ()9.2.5cp : To prove (|9.2.5ap . we first estimate ||92goo|k«- 
Using equation (|5.4.3ap . we have that 



(9.2.43) 

9*^500 - (.9°°)"'{ - S^'^aatgoo - agO^a^atgoo + Sffa^goo + 6H'{goo + 1) + Apo}. 

Using (|9.2.43p . Corollary EH Proposition EH the definition (|BT31) of Sat, 
Sobolev embedding, (|9.1.2ep . (|9.1.2hp . (|9.1.2dp . and (j9.2.4ap . it follows that 



(9.2.44) ||a2goo||ff«-i < Ce-'^^Sjv. 

We now use Corollary IA-31 and Proposition IA-6[ together with the definition 
of Sat, Sobolev embedding, (j9.1.2ap . (|9.1.2fp . (|9.1.2hp . fj.2.44p to obtain 

(9.2.45) 

||Pg,9a](ff00 + l)||L2 < ||.9™9a (^t'soo) " ^5 (5™a,^.9oo) |k« 

+ ||.9"'a5(9aafcffoo) - 95(5"'5,96500)lk« 

< CW^"\\H'^-i\\d^goo\\H--^+C\\d3'''\\H--A\dadbgoo\\H--^ 
+ C\\dg'"'\\H--4dadtgoo\\H--^ 
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This completes the proof of (|9.2.5ap . Inequahties (|9.2.5b[) and (|9.2.5cp can be 
proved similarly; we omit the details. 

□ 



10. The Equivalence of Sobolev and Energy Norms 

Our global existence proof is based on a standard strategy: showing that the 
Sobolev norms of Section [5] satisfy suitable bounds (they happen to be uniformly 
bounded for t > by Ce in the problem studied here), and then appealing to the 
continuation principle of Theorem l5.2l However, we do not have direct control over 
the growth of the norms; we can only control the norms indirectly through the use 
of the energies. In this section, we bridge the gap between the energies and the 
norms. More specifically, in the following proposition, we prove that under suitable 
bootstrap assumptions, the Sobolev-type norms and energies defined in Section [6] 
are equivalent. 

Proposition 10.0.2 (Equivalence of Sobolev Norms and Energy Norms). 

Let N > 3 be an integer, and assume that the bootstrap assumptions (|7.0.3ap - 
(j7.0.3cp hold on the spacetime slab [0,T) x for some constant ci > 1 and for 
11 ~ ^min- Let (6, y) be any of the pairs of constants given in Definition \6.2.1[ and 
let C(y) be the corresponding constant from Lemma \6.2.1\ There exist constants 
e'" > and C > depending on N, ci, r\min, Y, o,nd 6, such that i/ Sat < e'", then 
the following inequalities hold on the interval [0, T) for the norms and energies 
defined in (|6.1.1ap - (|BX^ . (|6.2.8ap - (|6.2.8dp . (|6.3.15p . and (PXT|) ; 

(10.0.46a) 

C-W\dtv\\L- + e-''\\dv\\L- + CmIHl-^) < Siy.8)[v,dv] < C{\\dtv\\L- + e-''\\dv\\L'^ + CmMl^) 
(10.0.46b) C^^Eggg+i-N < Sg„„+i-N < CEggg+i-N, 

(10.0.46c) C-'Eg„.N < Sg„,.N < CEg,,.N, 

(10.0.46d) C-'Eh„,N < Sh„,N < CEh„,N, 

(10.0.46e) C-'Eg.N < Sg-N < CEg,N, 

(10.0.46f) C-^Ep^p,a';N < Sp-p^u';N < CEp 

(10.0.46g) C^^Eat < Sat < CE 



Proof. The inequalities in (|10.0.46ap follow from the definition (|6.2.2p of f (^^) [v, dv], 
(|6.2.4p . the definition (|6.1.3p of Sat, Sobolev embedding, and (|7.0.4dp . The in- 
equalities in (jl0.0.46bp - (jl0.0.46dp then follow from definitions (|6.1.1al) - (|6.1.1cp . 
definitions (|6.2.8ap - (|6.2.8cp . and (|10.0.46ap . 

To prove (jl0.0.46f|) . for notational convenience, we set ii^ dsu^ , P == ds{P—p), 
and as in (j6.3.5p . we define vP =' —^^UaU"'- We now recall the definition (j6.3.13p of 
the energy current component J*^ associated to P, ii^ : 

(10.0.47) i° J^^P' + ^u'P + ^-^^f^g.pn'^u^. 
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Using the fact that < < i, (|9.1.1ap . (|9.1.1bp . (|9.1.1cp . (|9.2.2cp . (|9.2.2cp . 
( |9.2.2gp , Sobolev embedding, and the algebraic estimate |2m°P| < jy^^^P"^ + 

■(^-^Uav!^^ , it foUows that 



(10.0.48) 

-■^ 2^ — MaM 1^ |500| — Wait 

2(1 + c2)PuO 1 



— goau'^Ubu'' 

where W [P ~p, u^, v? , u^)^ . Integrating inequality (|10.0.48p over T"^, it follows 
that 

(10.0.49) \\j"\\L^>c{\\P\\h+e'-''j2\\^''\\l^}- 

a=l 

Similarly, it can be shown that 

(10.0.50) <c{||P||i.+e2^^^||^''||i.}. 

a=l 

Summing over all derivatives 9qW with \d\ < N, we have that 

(10.0.51) c-i {ll5<3(^-p)lli^ + ^'''Ell5«"1li4 

|q|<JV a=l 

< \\j°[d3W,dsW]\\Li 

\S\<N 

\a\<N a=l 

Inequality (|10.0.46f|) now follows from (|10.0.5ip and the definitions of Ep^p^u'-N 

and Sp-p,u'-N- 

Finally, (|10.0.46ep and p0.0.46gP follow trivially from definitions (|6.1.1dp . 
(j6.2.8dp . and (|6.4.ip . and from the previous inequalities. 

□ 

Remark 10.0.2. The pointwise positivity (|10.0.48p of J° is no accident. See 
[Chr07bj and |Spe09a| for further discussion concerning the coerciveness properties 
of the current J^. 
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11. Global Existence and Future Causal Geodesic Completeness 

In this section, we use the estimates derived in Section |9] to prove our main 
theorem. More specifically, we show that the modified equations (|5.4.3a|) - (|5.4.3e|) 
have future causally geodesically complete solutions for initial data near that of 
the background FLRW solution (g^i,,P,u^) on [0, cxd) x T'^, which is described in 
Section [4] We emphasize that this aspect of the theorem concerns only the mod- 
ified equations, and does not necessarily produce a solution to the Euler-Einstein 
equations p.2.4|) + p.2.6a|) - (j3.2.6b[) . However, as described in Section [5T6| if the 
Einstein constraint equations and the wave coordinate condition = are both 
satisfied along the Cauchy hypersurface 'S, = {x € Ai \ t = 0}, and the data for the 
modified equations are constructed from the data for the unmodified Euler-Einstein 
equations as described in Section [5731 then the solution to the modified equations 
is also a solution to the Euler-Einstein equations. The main idea of the proof is to 
show that the energies satisfy a system of integral inequalities that forces them (via 
Gronwall's inequality) to remain uniformly small on the time interval of existence. 
By Proposition 110.0.21 the norms must also remain uniformly small. The continu- 
ation principle of Theorem 15.21 can then be applied to conclude that the solution 
exists globally in time. 



11.1. Integral inequalities for the energies. In this section, we derive the sys- 
tem of integral inequalities that was mentioned in the previous paragraph. 

Proposition 11.1.1 (Integral Inequalities). Assume that {g fj_i, , P, ) , {^^v = 
0,1,2,3), (j = 1,2,3), is a classical solution to the modified system ()5.4.3a|) - 
(|5.4.3ep on [0,r) x and that the bootstrap assumptions (j7.0.3ap - (|7.0.3cp hold on 
[0, T) xT"^. Then there exists a constant e"" with < e"" < 1 and a uniform constant 
C depending only on N such that if S^it) < e"" holds on [0,2") and ti G [0,r), 
then the following system of integral inequalities is satisfied for t G [ti^T) by the 
energies defined in Section\Bi 



(11.1.1a) 



Ul_^{t) < C/^_i(ii) 



<o 



+ 




2[k-1 + q) e(i+'')"(^)C/2,_i(r) + CE^(t)[/^_i (r) dr, 



(11.1.1b) 
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(11.1.1c) 



(11.1. Id) 



I -4qHEl^.j,{T) + CEh,,,N{T)Eg„,,M{T) + Ce-'i"^^N{T)Eg„,,M{T 
Jti 



(ll.l.le) 



El,.Mi) < El,-M*i) + I He-'^"^El^.^^{r) + Ce~'^"^-E^{r)E^„.Mr)dT. 

Proof. To prove pi.l.lap . wc simply use (|6.3.16p . (|9.2.8f|) . and Proposition 110.0.21 
to conclude that 



(11.1.2) 



<0 3 



I {U%_,) < 2 (>. - 1 + g) a;e(i+')^^C/2,_, + 2e^^+''^''U^^, ^ II ||^«-. 

a=l 

<0 



< 2 - 1 + g) we(i+9)^^C/^_i + CEjvf/jv-i. 



Inequality (|ll.l.lap now follows from integrating (|11.1.2p from ti to t. 

To prove (jll.l.lbp . we make use of the following inequality, which we prove 
below: 



(11.1.3) ||a^(>[asW,asW])||^, <Ce-«"E2,. 

We also recall inequality (|6.3.17p : 



\a\<N 



Inequality ([ll.l.lbp follows from integrating (|11.1.4p from ti to t and using (|11.1.3p . 

To prove (|11.1.3p . we first recall equation (|6.3.14p . where W = (P, ii^^v?, v?)'^ 
dsW = {ds{P — p), dsu^ , dsu'^ , dsu^)^ , and as in (|6.3.5p . m° —-^Ua'u,"' ■ 
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(11.1.5) 



dt 
dt 

da 



(1+C2)P 

(1 + c- 



da 



p2 



(1 + c2)Pm° 



(1 



-{dagmW) 



(1+C2)P 

5oo(^°)2 + 2.goaM°u" 
(l + c2)Pu° 



gabU'^u'' 



9abU"'u'' 



-b 



(1 + c2)Fw'' 



{dagim)u''u" 



(1 



{dagob)u ii' 



(5t5oo)(^t°)' 



2(1 + c2)PmO 



Mo 

,2' 



2{1 + ci)P 



{dtgoa)u''u'' -r - 
(l + c2)P 



(l+c2)P(u"-l) 



{dtgab)u'^u 



:.b 



(>< - l)ujgabu''u + 



(9t5afc - 2ujgab)u''u^ 
2^ ^ 



(1 + C2)P 



<0 



2(1 + c2)P 



o„-,o 



4(1 + c2)P 



In the above expression, the terms ~ and © = are given by (|6.3.11bp and 
(|6.3.12p . The estimate (|11.1.3p now follows from Proposition 110.0.21 and (|11.1.5p 
using the following four steps: (i) we drop the negative (since k = 3c^ < 1 and 
gab is positive definite) fifth- from-last term on the right-hand side of ([ll.l.Sp : (ii) 
we bound each variation u*^, ^s, and ©^ in L^, and make use of (|9.2.9p . (|9.2.10ap . 
and (|9.2.10bP : (hi) we bound all of the remaining terms in L°°, and use Sobolev 
embedding together with the estimates (|9.1.3bp . (|9.1.3cp . (|9.2.2cp . (|9.2.2ep . (|9.2.3bp . 
(|9.2.3dp . (|9.2.3f|) . (|9.2.3hp . and ( |9.2.3jP ; (iv) we make repeated use of the estimate 
||uit'2'U3||Li ^ Iji'i ||i^2||l2 ||ti3||x,2, whcrc V2 and are terms estimated in step 
ii), and vi is estimated in step iii). 

To p rove (jll.l.ldp . we apply Corollary |6T3l using (I9.1.6ap . (I9.2.4bp . (j9.2.5bp . 
(|9.2.6bp , and Proposition 110.0.21 to estimate the terms on the right-hand side of 
(|6.2.11bp . and using definition (|8.1.2p to deduce that 2((7— 1) — t)o* < —4(7, thereby 
arriving at the following inequality: 



(11.1.6) jX^I,'n) < -MHEl^.^^ + CEH„;NEg,,,,N + Ce-^^^^NEg^^.N. 

Inequality (jll.l.ldp now follows from integrating (|11.1.6p in time. Inequalities 
(jll.l.lcp and (lll.l.lep can be proved similarly; we omit the details. □ 

Remark 11.1.1. The term CEh^^^.^Egg^.i^ in inequalities ([ll.l.ldp and (jll.l.6p 
arises from the Sgg^-Nj2'^=i^'''^~^''^^\\9"'''^a.jb\\H''^ term on the right-hand side of 
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(|6.2.11b[) . This term is dangerous in the sense tliat it docs not contain an exponen- 
tiaUy decaying factor, and looks like it could lead to the growth of Eg^^-N. However, 
as we shall see in the proof of Theorem 111.11 there is a partial decoupling in the 
integral inequalities in the sense that the CEh._^;N factor in the dangerous term can 
be controlled independently using inequality (jll.l.lcp alone. We will then insert 
this information into inequality (|ll.l.ld|) . and also make use of the negative term 
—AqHEg^ .jy to obtain a bound for Eg^^.]^. 

For completeness, we state the following version of Gronwall's inequality; we 
omit the simple proof. We will use it in Section [11.21 

Lemma 11.1.2. Let b{t) > be a continuous function on the interval [ti,T], and 
let B[t) be an anti- derivative ofb{t). Suppose that A > and that y[t) > is a 
continuous function satisfying the inequality 

(11.1.7) y{t)<A+ f b{T)y{T)dT 

Jti 

for t G [ii, r]. Then for t G [ti, T], we have 



(11.1.8) 



y(t) < Aexp Bit)-B{ti] 



□ 



In addition, in Section 111.21 we will apply the following integral inequality to 
(jll.l.ldp in order to estimate the energy Egg^.N{t). 

Lemma 11.1.3. |RS09| Lemma 11.1.3] Letb{t) > be a continuous non- decreasing 

function on the interval [0,r], and let e > 0. Suppose that for each ti G [0,T], 
y{t) > is a continuous function satisfying the inequality 



(11.1.9) 



yHt)<y^(t,) + 



-b{T)y^{T) + ey{T) dr 



for t G [ti,T]. Then for any ti,t Cz [0,T] with ti < t, we have that 



(ll-i-iO) y{t)<y{ti) + j^. 

b{ti) 

□ 

Proof. Let C be the "highest" curve in the (t, y) plane on which the integrand 
in (|11.1.9P vanishes; i.e. C = {{t,y)\y = ^^^}. Then by (|11.1.9p . above C (i.e. 
for larger y values), y(t) is strictly decreasing. Let y{t) achieve its maximum at 
tmax e [ti,T]. We separate the proof of (jll.l.lOp into two cases. Case i) assume 
that t,nax = tl. Then y{t) < y{t,nax) = y{ti) for t G [ti,T], which implies (|11.1.10p . 
Case ii) assume that t^ax G (^i, T]. Wc claim that y{tmax) < yj^ — j- For otherwise, 
the point [tmax, y {tmax)) Uss abovc C. Since y{t) is then strictly decreasing in a 
neighborhood of tmaxi it follows that there are times t* < tmax, with G (ii,T), 
at which y(t^,) < y(tmax)- This contradicts the definition of tmax- Using also the 
fact that is non-increasing, it follows that y{t) < y{tmax) ^ f^jj^ — y ^ fc(f7)' ^^^^ 
concludes the proof of pi.l.lOp . □ 
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11.2. The global existence theorem, hi this section, we state and prove our 
main theorem, which provides global existence criteria for the modified equations 
()5.4.3a|) - (|5.4.3ep . and future-stability criteria for the unmodified Euler-Einstein 
equations (|3T4l) + (|3.2.6ap - ()3.2.6b|) . 

Theorem 11.1 (Future-Stability of the FLRW Family). Assume that < 
< 1/3 and N > 3. Let g^^ = gfj.iy\t=o, 2/<'^^ = dtg^^\t=o, (Mi = 0,1,2,3), 
P = P\t=o = p\t=o, = u^t=o, (j = 1,2,3), gaf3U°'u'^\t=o = -1, be initial data 
(not necessarily satisfying the wave coordinate condition or the Einstein constraints) 
on the manifold T'^ for the modified Euler-Einstein system ()5.4.3a|) - ()5.4.3e|) . and 

let Sn = Sg._N + Sp-p_u';N + Un~i be the norm defined in ()6.1.3p . Assume that 
there is a constant ci > 2 such that 

(11.2.1) -SabX^'X'' KgabX'^X'' K^SabX'^X'', \/{X\X^X'') eM.\ 
Cl 2 

Then there exist a small constant eo with < eo < 1 o,JT-d a large constant 

C* such that i/ e < eo "■iT'd Sn{0) < C^^t, then the classical solution [g^y,P = 

g3(i+Cs)f2p^ provided by Thcorem \5.1\ exists on A4 == [0, oo) x T'^, and 

(11.2.2) SNit)<e 

holds for all t > 0. Furthermore, the time Tmax from the hypotheses of Theorem \5.S\ 
is infinite, and the spacetime-with-boundary {M,gfj,u) is future causally geodesically 
complete. 

Finally, if the initial data {i/f^i,, K^^, P — p,u^) for the modified system are 
constructed from initial data {gj^,K_ji^,p,]M'), {j,k ^ 1,2,3), for the unmodified 
Euler-Einstein equations p.2.4p + p.2.6ap - (j3.2.6bp as described in Section \5.3\ 
then the solution {g^u^P = e~^'^^^'^'^^P,u^) to the modified system is also a future 
causally geodesically complete solution to the unmodified equations. 

Remark 11.2.1. It is possible to restate the stability criteria in terms of quantities 
that manifestly depend only on the closeness of the initial data (ir^,5^.j,,^jfe,P,M-') 
for the unmodified system to the corresponding data for the FLRW background 
solution (g,p,u). For example, a sufficient condition for global existence and future 
causal geodesic completeness would be 

3 3 

(11.2.3) iii.fe ~ ^.-^11^^"+^ + E wkjk-^mkU- 

3 

+ \\p~p\\h^ + Y,\\^'\\h'' <e, 

where e is sufficiently small. This is because the condition (|11.2.3p implies that 
S7v(0) < Ce and furthermore (by Sobolev embedding) that a condition of the form 
()11.2.ip holds (i.e. that the hypotheses of Theorem UTTTl hold) . 

To see that SAr(O) < Ce follows from (|11.2.3p . wc first use the definition (|6.1.3p 
of Sn (0) , the construction of the modified data described in Section 15.31 and the 
triangle inequality to deduce that 
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(11.2.4) 

3 3 

s^(o) < 2\\M^)-r^k^,\\HN + j2 mjH-+2 \\^m^,-K,k\\H- 

3 3 

< 2||(f " - 5'^')k^,\\H« + 2W\k^, - c.(0)<5,,) 

3 3 

+ E + 2^(0) E Wgjk^^AH- 

j,k=l j.k=l 

3 3 

j,k=l 3 = 1 

3 

Now using Corollary I A-4[ Proposition IA-5[ and Sobolev embedding, it follows that 
if (|11.2.3I) holds and if e is sufficiently small, then the right-hand side of (111.2. 4p is 
< Ce. 

Proof. See Remark 19.2.11 for some conventions that we use throughout this proof. 
We only discuss the issue of global existence and obtaining the uniform bound 
SAr(t) < e for t S [0,00). The future causal geodesic completeness of the result- 
ing spacetime-with-boundary follows from this bound via the arguments given in 
jRin08| Propositions 3 and 4] and jRS09| Theorem 11.2]. Our global existence proof 
relies upon a standard bootstrap-style argument that ultimately relies on Theorem 
15.21 i.e., we will make assumptions concerning the size of the energies and concern- 
ing (7^1,, and we will use these assumptions, together with assumptions on the data, 
to deduce an improvement. In effect, we will avoid the four breakdown possibilities 
of Theorem 15.21 which will allow us to conclude global existence. 

To begin the analysis, we invoke Theorem 15.11 which shows that if e is small 
enough and SAr(O) < e, then there is a local solution (5^1,, existing on a 

(non-trivial) maximal interval [0, T) on which the following bootstrap assumptions 
hold: 



(11.2.5) Sjv(i)<e, 

(11.2.6) |goO + 1| < Tlmiri, 

(11.2.7) c^HabX'^X'' < e-^''gabX''X'' < c^5abX''X\ y{X\X^, X^) e M.\ 

3 

(11.2.8) Elffoal' <Tl™:nCr^e2(i-«)^\ 

a=l 

(11.2.9) F>0. 
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Observe that the rough bootstrap assumptions (|7.0.3ap - (|7.0.3cp are included in 
the above assumptions. By maximal interval, we mean that 

(11.2.10) 

T = sup {i > I The solution exists on [0, i] x T^ and (|11.2.5|) - (|11.2.9|) hold}. 

We may assume that T < oo, since otherwise the theorem follows. The remainder 
of this proof is dedicated to reaching a contradiction if e is small enough and C* is 
large enough. 

Our first assumption, which we use repeatedly throughout this proof, is that e 
is small enough so that Propositions I9.1.T1 19.2.11 and I10.0T21 are valid on [0,T). 
We will make repeated use of Proposition 110.0.21 throughout this proof without 
explicitly mentioning it each time. 

Wc now address the bootstrap assumptions (|11.2.6|) - (|11.2.9|) , with the intent of 
showing an improvement. First, we note that the assumption Sat < e implies that 

(11.2.11) ||at(e-2%.fc)|li^ = Wdthjkh^ < Cee-'^"K 

We then use (|11.2.1ip to integrate in time from t = Q, concluding that 

(11.2.12) ||e-2ng^.,(i,.)-5^.,(.)||^^ <Ce. 

By (|11.2.ip and (|11.2.12p . it follows that if e is small enough, then on [0,T) x T^, 
we have that 

(11.2.13) —SabX^X'' < e-'^'^gabX^X'' < ^dabX''X\ ViX\X^, X^) G 

Furthermore, by the definition (|6.1.3p of Sw, Sobolev embedding, and (|11.2.5|) . the 
following inequalities hold on [0,T) : 

(11.2.14) ||goo + l||L- <Ce, 

(11.2.15) |lffo,|U== <Cee(i~^)^*, 

(11.2.16) IIP-pIIl- <Cee-«^*. 

Inequalities ()11. 2.131 - (|11.2.16|) show that if e is small enough, then the bootstrap 
assumptions (|11.2.6|) - (|11.2.9|) can be strictly improved on the interval [0,T). 

To complete our proof of the theorem, we will show that if e is small enough and 
C* is large enough, then the bootstrap assumption (|11.2.5p can be improved by 
replacing e with e/2; the primary tool for deducing an improvement is Proposition 

111. 1.11 Throughout the remainder of the proof, we set SAr(O) = e. To be gm our 
proof of an improvement of (|11.2.5|) . we use a very non-optimal application of 
Proposition 1 1 1 . 1 .T] with ti = 0, deducing (with the help of Proposition [TP. 0.2[) that 
on [0,T), wc have that 

(11.2.17) S%it)<CS%iO)+ f cS%ir)dT. 

Jt=0 
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Applying Lemma [11.1. 21 (Gronwall's inequality) to (|11.2.17p . using Sjv(O) = e, and 
using Proposition 110.0.^ we conclude that the following weak inequalities hold on 

[o,r): 

(11.2.18) Sjv(t) < Cke"', 

(11.2.19) EAr(i) < Cee=*. 

Remark 11.2.2. The weak inequality pi.2.18|) already shows that the time of 
existence is at least of order c~^|ln(Ce)|, if e is sufficiently small. 

We now fix a time ti G [0, T); ti will be adjusted at the end of the proof. Roughly 
speaking, it will play the role of a time that is large enough so that the exponentially 
damped terms on the right-hand sides of the inequalities of Proposition 111. 1.11 are 
of size ^ e. 

We now estimate the energy Eggg^i-N{t). To begin, we drop the negative term 
—AqH E'^^^j^^.pj{t) on the right-hand side of inequality (jll.l.lcp and use the boot- 
strap assumption (|11.2.5p to deduce that 

(11.2.20) El^+iA^) < Ko+iMt^) + /* e-^^^ dr. 

Using (|11.2.19| ) at time ti and performing the integration on the right-hand side of 
(jll.2.20|) . we have the following inequality for t S [ti,T) : 

(11.2.21) Eg,,o+i;N{t) < C{le^*' + ee-«^*i/2}. 

Also using (|11.2.19|) to obtain an upper bound for Eggg^i-]^{t) on [0,ti], we have 
established the following inequality, which is valid for t S [0, T) : 

(11.2.22) i?soo+i;A^W < C{le''' + ee-""''/^}. 

To estimate Eh,,-N, we can appeal to inequality ()ll.l.le[ ) argue as we did for 
£"900+1:^, thus obtaining that the following inequality holds for t G [0,T) : 

(11.2.23) Eh„-Nit) < C{lef'' + ee-'^"'''^). 

To estimate Egg^-]^{t), we use (jll.l.ldp . the bootstrap assumption ()11.2.5p . and 
(|11.2.23p to arrive at the following inequality valid for t G [ti,T): 

(11.2.24) 

El,,At) < El,,M + / -MHEl^.^^ir) + C{ee^*^ + ee~^"'-'^}Eg,,,N{r) dr. 

JT = ti 

Applying Lemma [TTX31 to (|11.2.24p . with y{t) = Egg^.^it) and b{t) = AqH in the 
lemma, and also using (jll.2.19p at time ti, we conclude that the following inequality 
holds on [0,r) : 



(11.2.25) 



Ego,;N{t)<C{ee'^'^+ee-'^"'^/^}. 



Stability of the FLRW Family 



69 



To estimate UN-i{t) on [ti, T), we first use (|ll.l.lap and (|11.2.5p to deduce that 

for t G [ti,T), we have 

(11.2.26) 

t ^° 

C/^-i(i) < f/^-i(ii) + / 2(^^-1 + g)e(i+?)^-[/2,_^(r) + Cet/jv-i(T) dr. 

Applying Lemma [TTX^ to (|11.2.26p . with y(t) = ?7jv-i(i) and fe(t) = |2(>< - 1 + 
(3')|e^^+'?-'^* in the lemma, and also using ()11.2.19p . it follows that on [0, T), we have 

(11.2.27) UM-i{t) < C{ee'=*i + ee-'^"''/^}. 

To estimate Ep-p^u';N{t) on [ii,T), we simply use (|11.2.5p and (|11.2.19p to 
estimate the two terms on the right-hand side of (jll.l.lbp : 

(11.2.28) 

Also using (|11.2.19p to estimate Ep-p^u';N{t) on [0,ti], we conclude that the fol- 
lowing inequality is valid on [0, T) : 

(11.2.29) Ep^p,^.,N{t) < C{le'*' + eer""''/^}. 

Adding (|11.2.22p . (|11.2.23p . (|11.2.25p . (|11.2.27p . and (|11.2.29p . referring to defi- 
nition (|6.1.3p . and using Proposition 1 10. 0.21 it follows that on [0,r) 

(11.2.30) SAr(i) < C{ee='i + ee-«^*i/2}. 

We now choose <i such that Ce"'^*^/^ < -j, and I such that eCe'^^^ < je, where 
the constant C is from the right-hand side of (|11.2.30p . This implies (with the help 
of Proposition 1 1 . ."2| that on [0,r), we have that 

(11.2.31) SN{t) < ie. 

We remark that in order to guarantee that the solution exists long enough (i.e. that 
T is large enough) so that ti G [0, T), we may have to further shrink e; see Remark 
111.2.21 We also remark that by the above reasoning, it follows that the constant C* 
from the conclusions of the theorem can be chosen to be ACe'^^^ , where C is from 
the right-hand side of (|11.2.30p . 

Combining (|11.2.13p - (|11.2.16p and (|11.2.3ip . and using Sobolev embedding, 
it now follows that none of the four existence-breakdown scenarios stated in the 
conclusions of Theorem 15.21 occur. Using the continuity of Sjv(i), it also follows 
that the solution can be extended to an interval [0,T + 5] on which the bootstrap 
assumptions pi.2.5p - (|11.2.9p hold. This contradicts the maximality of T and 
completes the proof of the theorem. 

□ 
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12. ASYMPTOTICS 

In this section, wc provide a theorem that strengthens the conclusions of Theorem 
111.11 More specifically, we show that g^ji,, g^"^, P ~ p, u^, (/i, — 0,1,2,3), and 
various coordinate derivatives of these quantities converge as t —> cxd. Furthermore, 
some of the decay rates (for example, those in (|12.0.34a|) and (jl2.0.34bp ) are a 
significant improvement compared to the rates that can be directly obtained from 
the bound Sat < e, which was derived in Theorem lll.il The results of this section 
parallel the ones obtained in |Rin08| Proposition 2] and |RS091 Theorem 12.1]. We 
remark that they are not optimal, and more information could be extracted with 
additional work. 

Theorem 12.1 (Asymptotics). Assume that the initial data (^^i/, i^^i/, P, w^), 
{fj.ji' = 0,1,2,3), {j = 1,2,3), .ga^M"u''|t=o = —1, for the modified Euler- Einstein 
system (|5.4.3ap - (|5.4.3ep satisfy the assumptions of Theorem \11.1\ including the 
smallness assumption SAr(O) < C~^e, where < e < ep. Let g^^^ denote the inverse 
of g^iu- Assume in addition that N > 5, and let {g^^, P,u'^) be the future-global 
solution launched by the data. Then there exists a constant ti satisfying < £2 < co 
such that if e < £2, then there exists a Riemann metric = 1,2,3J, with 

corresponding Christoffel symbols r|°^'', {i,j,k = 1,2,3), and inverse <7(^) on T'^, 
a function P(oo) on T^, and a time-independent vectorfield ^^l^-) on such that 

- 93k e H^, g[l^ - g^^ e i?^, - p £ H^-\ u\^^ e and such 

that the following estimates hold for all t > : 



(12.0.32a) 
(12.0.32b) 



llff(l) 



g^^Wn- < Ce, 



(12.0.33a) 
(12.0.33b) 
(12.0.33c) 
(12.0.33d) 
(12.0.33e) 
(12.0.33f) 
(12.0.33g) 
(12.0.33h) 



-2n 



'93>^-g^T^\\H^<Cee-^"' 



lie gjk-g'kl 



V''^-.C)ll 



glL^\\H«<Cee-'^"' 



''dtg.k 



We^^'dtg^" - 



,ll^U.-.<Cee-'^^ 
-2..4r^||^«<Cee-'^* 
2ug<^f\\H^-.<Cee-'"' 



(12.0.34a) 
(12.0.34b) 
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(12.0.35a) 
(12.0.35b) 
(12.0.35c) 
(12.0. 35d) 

(12.0.36a) 
(12.0.36b) 



11.900 + 1||h" <Cee-«^*, 
i3oo + 1||h«-2 <Ce(l+t)e-2ff*, 



\\dtgm + 2w(goo + l)||if«-2 < Cee~^"\ 



-2n 7 



I 



jk 



II 



1 < Cee-«^*, 
< C(l +t)ee-2-f^*. 



In the above inequalities, Kjk is the second fundamental form of the hypersurface 
{t = const}. 



Furthermore, 


we have that 








(12.0.37a) 






M\h- 




(12.0.37b) 








< Cee-2(i— )^^*^ 


(12.0.37c) 




i 

(oo) 






(12.0.37d) 








< Ce, 


(12.0. 37e) 


11^ 






< Cee-""', 


(12.0. 37f) 


11^ 


- ^'(oo) 






(12.0.37g) 


ll^^( 


oo) 




< Ce. 



Proof. Wc only provide a sketch of the proof, since most of the details can be 
found in the proof of |RS091 Theorem 12.1]. The main idea behind the improved 
decay rates is that many of terms in the modified equations (|5.4.3ap - (|5.4.3eP can 
be treated as inhomogeneities for equations that have a more favorable structure. 
In fact, many of the terms that are of principal order from the point of view of 
the number of derivatives can be treated as lower-order terms in the sense of decay 
rates. Of course, this approach is only viable for treating the lower-order derivatives 
of the solution. A related consequence is that some of our estimates are proved only 
in Sobolev spaces of lower order than the ones the data belong to. 

The estimates (|12.0.32ap - (|12.0.36bj) can be proved by a straightforward adap- 
tation of the proof of |RS09| Theo rem 12.1], and we omit the details. Some very 
important ingredients in the proof of jRS09| Theorem 12.1] were the following up- 
graded estimates 



(12.0.38) 


Wdtgoj 




<Ce, 


(12.0.39) 






< Ce, 


(12.0.40) 






< Cee", 


(12.0.41) 






< Cee-^^ 


(12.0.42) 


\\dthjk 


\hn-2 




which we will use below. 









To prove (|12.0.37cp . we first use (|5.7.5cP to deduce that 
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(12.0.43) (9t(e(2— ^ e'-^-'^^^A'^. 

Wc now introduce the following non- negative fluid energy 'au*-N-2{t), which is de- 
fined by 

(12.0.44) €^;N-2'^j2\\e^'-''^''u^rH^-^- 

Wc then differentiate under the spatial integrals implicit in the definition p2.0.44p 
and use (|12.0.43p plus the Cauchy-Schwarz inequality for integrals, thus arriving at 
the following inequality: 

(12.0.45) |(e;A'-2) <2<§'„.;^-2e(^-)^lA'-'||^«-.. 

Using Corollarv lA-31 Proposition lA-5[ the definition (|6.1.3p of Sat (which satisfies 
< e), the definition (|12.0.44p of S'u'-n-2, Sobolev embedding, the estimates of 
Proposition 19. 1.1 1 and Proposition 19.2.11 (|12.0.33ap - (|12.0.35dp . and the improved 
estimates (|12.0.38p - (|12.o!42p . it follows that 

(12.0.46) ||A'J'||h«-2 < Cee-^^ + C<^u':N-2e^^''~^'>^ ■ 

We remark that the CS'u';N-2e'^^'^~'^^" term on the right-hand side of (|12.0.46p 

2 

arises from the jj^^^^jpH'^^ A' term on the right-hand side of (|5.7.6cp : this term 
contains terms that are the product of the spatial metric components gat with the 
cube of the spatial four- velocity components . We now insert the bound p2.0.46p 
into the right-hand side of (|12.0.45p . use the initial condition /^u' -n -2{0) < Ce, and 
apply Gronwall's inequality, thus concluding that 

(12.0.47) ^„.;w-2(t) < Ce. 

Inserting the bound (|12.0.47p into the right-hand side of (jl2.0.46p and revisiting 
(|12.0.43p . we have that 



(12.0.48) ||5t(e(2— )^^u^')||jj„_2 < Cee-''"' + Cee^^''-^'>"K 



From (jl2.0.48|) . it easily follows that there exist functions ul Ax^^x^.x^) such that 



(12.0.49) ||e(2— - u]. \\hn-2 < Cee""^* + Cee^^^"!)^*, 



(oo) 

where 



(12.0.50) \\ul^^\\H--2<Ce. 

We have thus shown (|12.0.37cp and (|12.0.37d|) . 

The estimates (jl2.0.37ap - (|12.0.37bp and (|12.0.37cp - ( |12.0.37gP can be proved 
similarly with the help of equations (|5.4.4I) and (|5.7.5ap : we omit the details. 

□ 
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Appendices 

A. SOBOLEV-MOSER INEQUALITIES 

In this Appendix, we provide some standard Sobolev-Moser estimates that play a 
fundamental role in our analysis of the nonlinear terms in our equations. The propo- 
sitions and corollaries stated below can be proved using methods similar to those 
used in |II6r971 Chapter 6] and in jKM81| . The proofs given in the literature are 
commonly based on a version of the Gagliardo-Nirenberg inequality |Nir59| , which 
we state as Lemma FA- li together with repeated use of Holder's inequality and/or 
Sobolev embedding. Throughout this appendix, we abbreviate ~ LP{T^), and 

Lemma A-1. If A'I,N are integers such that < M < and v is a function on 
such that v e Hl^ < oo, then 

(A.l) \\d^''^v\\L..,M < C(M, 7V)||H|^#||9W,;|||. 

Proposition A-2. Let M > be an integer. If {va}i<a<i are functions such that 
Va S ||9^^^''vQ||i2 < oo for 1 < a < and (5?i,--- , (5; are spatial derivative 
multi-indices with \di \ + ■ ■ ■ + \di \ = M, then 

I 

(A.2) \\{ds,v,){d3,V2) ■ ■ • {ds,vi)\\L- < C{1, M) {\\d^''^Va\\L^- n WvbWL^). 

a=l b^a 

Corollary A-3. Let M > 1 be an integer, let ^ be a compact set, and let F G 
C*^(.S) be a function. Assume that v is a function such that v{T^) C K and dv G 
H^-\ Then d{F o v) e H"'-^, and 

M 

(A.3) \\d{Fov)\\HM-^ < CiM)\\dv\\m'-^J2\^^'^\-^ML-- 

1=1 

Corollary A-4. Let M > I be an integer, let R be a compact, convex set, and let 
F G C^^ (M.) be a function. Assume that v is a function such that v{T^) C K and 
w — w G H^'^ , where v G R is a constant. Then Fov — Fov(^ H^^ , and 

(A.4) 

M 

\\Fov ^FovWhm < C{M){\F^'^U\\v - vh. + \\dv\\H^,-. ^ \F^'^U\\v\f^2}. 

1=1 

Proposition A-5. Let AI > 1,1 > 2 be integers. Suppose that {va}i<a<i o.i'^ 
functions such that Va G L°° for 1 < a < I, that vi G H'^^ , and that di'a G II^^~^ 
for 1 < a < I — 1. Then 
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(A.5) 

i-i i-i 

\\viV2---vi\\hm < C(/,M)|||z;,||^A. IKh^ + ^\\dva\\H'"-^ n ll^fclli^}- 

a— 1 a— 1 b^a 

Remark A.l. The significance of this proposition is that only one of the functions, 
namely w;, is estimated in L^. 

Proposition A-6. Let M > 1 be an integer, let ^ be a compact, convex set, and let 
F G C^^(i?) he a function. Assume that vi is a function such that vi{T^) C that 
dvi e L°°, and that d^^'^'>vi € L'^. Assume that V2 e that d^^'^^\2 S L^, and 
let d be a spatial derivative multi-index with with |d| = M. Then ds {{F o Vi)v2) — 
{F o vi)dsV2 G L'^, and 

\\ds{,{F ovi)v2) - {F o vi)dsiV2\\L-^ 
(A.6) 

M 

<C{M)[\F^^-^U\\&v,\M\d^^^^^^^ 

1=1 

Remark A. 2. The significance of this proposition is that the Af*'* order derivatives 
of V2 do not appear on the right-hand side of (jA.Sp . 
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